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1 Introduction generale 



La recherche dans le domaine des sciences des materiaux (developpement de la fabri- 
cation de structures artificielles en couches d'alhage . . . ) et en particuher dans celui 
des couches minces et multicouches magnetiques a pris une ampleur considerable 
grace a I'enjeu economique de I'enregistrement magnetique et magneto-optique. Leur 
interet s'est rapidement elargi avec le develop emment des couches minces deposees 
sur des substrats mono-cristallins. Les disques durs, bandes magnetiques et disques 
optiques sont aujourd'hui les principaux supports utilises dans les ordinateurs pour 
stocker des quantites considerables de donnees|B El- Les systemes a disques durs 
magnetiques representent aujourd'hui pres des deux tiers de I'ensemble du marche des 
peripheriques de stockage. La surface occupee par un bit (un domaine magnetique 
d'aimantation planaire) a ete reduite d'un facteur cent mille en trente annees de 
recherches et d'innovations. 

Les capacites toujours croissantes des ordinateurs ont stimule des recherches numeriques 
de plus en plus ambitieuses portant sur des systemes de taille croissante. Ces simu- 
lations peuvent etre confrontees a des solutions analytiques permettant d'en preciser 
les limites de validite, ou de servir de guide dans I'exploration de comportements 
mal ou peu connus. On dispose ainsi de moyens permettant de verifier les hypothese 
theoriques, de les modifier en fonction des donnees ou encore de decouvrir de nou- 
veaux territoires. 

La recherche de materiaux nouveaux de hautes performances, la fabrication de 
materiaux composites, la course a la miniaturisation en microelectronique, donnent 
un role de plus en plus grand aux surfaces et interfaces. La fabrication de couches 
tres fines de super-reseaux a ete proposee dans les annees 1970 par L. Esaki et R. Tsu 
des laboratoires IBM. C'est ensuite Y. Cho et J.R. Arthur des laboratoires Bell qui 
ont developpe la technique d'epitaxie par jets moleculaires qui permet de controler et 
d'interrompre la croissance a moins d'une monocouche atomique pres. Des proprietes 
interessantes resultent de I'augmentation de la densite des interfaces et du caractere 
bidimensionnel de ces materiaux et de la regularite de la structure artificielle multi- 
couche, laquelle est determinee par la nature (topologique et chimique) de I'interface 
entre les differents metaux. 

L'existence d'une surface pent modifier les proprietes d'un materiau en question, 
et c'est a travers cette surface qu'il interagit avec I'exterieur. D'un point de vue 
industriel, les etudes ont ete principalement motivees par leurs applications dans 
les domaines de I'enregistrement magnetique perpendiculaire (supports de stockage 
haute densite de I'information), de la magnetoresistance (capteur et tete de lecture) 
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1 Introduction generale 



de la supraconductivite (depot sur des cavites resonantes) ou de la tribologie (couches 
de protection contre I'usure). 

Les deux principales techniques utilisees dans le domaine de la simulation en 
physique statistique de la matiere condensee sont la Dynamique Moleculaire (DM) 
et la methode Monte-Carlo (MC). Ces simulations reposent sur une description 
moleculaire en termes de positions et quantites de mouvement d'un ensemble de 
particules qui constituent le systeme. L'approximation de Born-Oppenheimer est 
adoptee, c'est-a-dire que I'on exprime I'Hamiltonien du systeme en fonction des vari- 
ables nucleaires, le mouvement (rapide) des electrons ayant ete moyenne. Faisant 
l'approximation supplementaire qu'une description par la mecanique classique est 
adaptee, on pent ecrire I'Hamiltonien H du systeme comme la somme de I'energie 
cinetique K, et de I'energie potentielle U de I'ensemble des atomes en fonction des 
coordonnees r, et des quantites de mouvement pj de chaque atome i. 

H( r^, p^)= K( p^)+ U( r^) 

L'energie cinetique a une expression bien etablie alors que I'energie potentielle qui 
decrit les interactions inter-atomiques est en general mal connue et constitue, de ce 
fait, une limitation importante de ces methodes de simulation. Les potentiels inter- 
atomiques utilises souffrent souvent de la faon tres approchee et indirecte avec laquelle 
les electrons sont traites. A priori, I'effet des electrons est cache dans les parametres 
du potentiel, qui doivent etre determines par ajustement. En general, cet ajustement 
se fait sur des donnees experimentales, qui ne sont pas toujours facilement accessibles 
ou disponibles. 

Parmi tons les alliages a base de cobalt deja etudies et qui out suscite un re- 
gain d'interet de la communaute scientifique depuis le developpement de la fabri- 
cation de structures artificielles en couches minces, le systeme CoPt est I'un des 
plus prometteurs car il est susceptible de presenter de fortes energies d'anisotropie 
magnetocristalline du fait de la stabilite de la phase anisotrope LIq et du couplage 
spin-orbite tres eleve [3 1^ [3] . 

Ce memoire est organise de la maniere suivante: 

Dans la premiere partie de ce travail, nous donnerons un aperQu sur les methodes 
de simulation dans les materiaux; nous rappellerons les proprietes pertinentes du 
modele en Liaisons Fortes [HI et nous en deduirons I'expression de la densite d'etat 
electronique locale en volume qui permet le calcul de I'energie de bande par site[71IS]. 

Dans le chapitre 2, nous presenterons les methodes de calcul des parametres du 
potentiel inter-atomique dans l'approximation au second moment (SMA)'S' de la 
densite d'etats, qui permettent de reproduire raisonnablement les energies de forma- 
tion des differentes phases (LIq, A1 et LI2) de I'alliage CoPt. Nous determinerons 
ensuite les diverses energies de col du systeme, ces energies jouant un role important 
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en simulation Monte-Carlo. 

Dans le troisieme chapitre, apres quelques generalites sur la theorie des transitions 
de phases et une description de la methode Monte-Carlo (MC) ainsi que des differents 
algorithmes de simulation, nous developperons notre modele MC en Liaisons Fortes 
dans le cadre de I'approximation au second moment (MC-SMA). Ce modele sera 
utilise pour evaluer la transition ordre-desordre LIq — > Al de I'alliage CoPt. 

Dans le chapitre 4, nous rappellerons quelques generalites sur les surfaces et les 
methodes de caracterisation structurale des films minces, en particulier dans I'alliage 
CoPt. Nous determinerons ensuite les energies d'adsorption du Co et du Pt sur la 
surface de I'alliage CoPt avec et sans relaxation. Et finalement nous presenterons 
une etude preliminaire sur la simulation de marches. 

Enfin, nous terminerons par un travail qui s'inscrit dans la continuite des etudes 
de cinetiques d'ordre a longue distance, dans le cadre d'un modele d'interaction de 
paires, entreprises initialement par Yaldram et Kentzinger 0, dans la phase B2 en 
collaboration avec le L.P.C.Q de Tizi-Ouzou. L'etude s'est poursuivie, dans le cadre 
de deux theses de magister (A.Kerrache et M.Hamidi ^J^J) dans les structures 
ordonnees sur reseau cubique a faces centrees et hexagonales compactes. Le but 
essentiel de ce chapitre est d'etudier I'influence de I'epaisseur du film et des surfaces 
sur la temperature critique et I'energie de migration atomique, et de compare aux 
resultats obtenus precedamment. 
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2 Dynamique Moleculaire et Modele en Liaisons 
Fortes 

2.1 Introduction 

La physique de la matiere condcnsee ainsi que la science des materiaux sont con- 
cernees fondamentalement par la comprehension et I'exploitation des proprietes des 
systemes d 'electrons et de noyaux atomiques interagissant entre eux. Ces interactions 
sont assez bien decrites par la mecanique quantique, et presque toutes les proprietes 
des solides peuvent etre ainsi etudiees en utilisant des outils de calculs convenables. 
Les calculs des proprietes des solides peuvent etre effectues a I'aide d'une grande 
variete de methodes, de semi-classiques aux approches purement quantiques. 

L' extraordinaire montee en puissance des ordinateurs, qui va de pair avec leur 
vulgarisation massive, ne pouvait manquer d'influencer (pour ne pas dire boulcvcrser) 
les methodes d'etude en sciences et techniques, notamment celles necessitant des 
calculs tres longs irrealisables auparavant. La physique a ainsi vu se developper una 
activite de simulation numeriquc que Ton pcut situcr entrc la theorie pure (c'est-a- 
dire analytiquc) et I'experience et qui sort souvcnt de relais entre les deux. Dans 
cette partie, nous allons voir que la simulation numerique pent paraitre soit comme 
une experience, soit comme une approche theorique, qu'elle peux etre toujours reliee 
a une etude analytique des structures mathematiques sur lesquelles elle s'appuie et 
egalement, si possible, a de vraies experiences. 
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2 Dynamique Moleculaire et Modele en Liaisons Fortes 



2.2 Simulation de Dynamique Moleculaire 

D'apres la mecanique classique, la connaissance de la position et de la vitesse d'une 
particule donnee et des forces qui s'appliquent sur elle, permet de predire toute sa 
trajcctoire future. Cette conception a ete defendue par Helmholtz qui considerait que 
le probleme de la science physique consistait a ramener les phenomenes naturcis a des 
forces d'attraction et de repulsion dont I'intensite depend uniquement de la distance. 
Ces memes hypotheses sont admises lorsqu'on utilise un potentiel d'interaction pour 
etudier un phenomene physique. 

Les efforts ont sur tout porte jusqu'a present sur le choix des fonctionnelles ser- 
vant a definir le potentiel: elles font intervenir un nombre plus ou moins important 
de parametres ajustables. L'avantage principal dc ccttc technique est sa simplicite 
puisque du choix du potentiel et de son ajustement decoule ensuite une certaine 
qualite predictive. 

La dynamique moleculaire est actuellement la methode de simulation la plus utilisee 
en science physique (en matiere condensee). Elle consiste a etudier la trajectoire des 
atomes en interaction, et elle est basee sur la resolution numerique des equations de 
la mecanique classique newtonienne. Elle prend en compte I'interaction de chaque 
atome avec tous les atomes voisins formant le systeme pour le calcul des forces. La 
description des potentiels d'interactions varie d'une methode a I'autre: d'approche 
"ab-initio" a approche semi-empirique. 

La description newtonienne d'un tel systeme est donnee par: 



mj =masse de I'atome i 

Fi =force agissant sur I'atome i 

fi =coordonnees cartesiennes de I'atome i 



L'expression analytique de la force qui s'applique a chaque atome s'ecrit dans le 
cas de potentiels a deux corps: 




m 



(2.1) 




(2.2) 



3 



fij = distance entre les atomes i et j . 
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2.2 Simulation de Dynamique Moleculaire 



Lcs expressions de V(rjj) peuvent etre assez variees (voir table 1.1). Les forces 
electrostatiques font I'objet d'un traitement special par somme d'Ewald. 



Lennard-Jones 
(gaz rares) 


Lennard-Jones generalise 
p>q, (metaux) 


Buckingham 
(atomes a coeur "mou") 


C^.((f )^'-2(f )6) 


^..((f F-2(f )") 


Aje'^ - ) 



Table 2.1: Quelques expressions tres usuelles de potentiels interatomiques a deux 
corps. Les coefficients dependent des types d'atomes consideres et on ap- 
plique souvent des regies de combinaison comme: ^ab = V ^aa^bb 6t: 

Le succes d'une simulation par dynamique molculaire repose notamment sur la 
rigueur de la procedure d'integration de la loi de Newton pour chaque particule i de 

masse m^ et de position r^. 

L'integration d'un tel systeme d'equations differentielles se fait en subdivisant la tra- 
jectoire en une serie d'etats discrets, separes par des intervalles de temps tres court, 
dont la duree definit le pas d'integration. 

Les methodes numeriques utilisees pour resoudre des systemes d'equations differentielles 
du type de I'equation (1.1) presentent une grande diversite; cependant toutes consis- 
tent en approximation de l'integration sur une variable continue par une sommation 
sur une variable qui prend des valeurs discretes du type: 

ti = {i- l)dt 

ou dt est le pas d'integration dont la valeur est petite mais finie. 

Cette integration est realisee par diverses methodes mathematiques. On peut citer 
la methode de Range-Kutta Gill, la methode des differences centrales, la methode 
d'Euler-Cauchy, methode de Verlet 

II faut choisir un algorithme dc resolution numerique qui verifie les lois fondamen- 
tales de conservation de I'energie, de I'impulsion totale et du moment angulaire. Get 
algorithme doit egalement etre assez rapide. Plusieurs algorithmes ont ete testes: 
I'algorithme d'Euler, I'algorithme de Runge-Kutta d'ordre 2 (RK2), I'algorithme de 
Runge-Kutta d'ordre 4 (RK4) et I'algorithme Runge-Kutta d'ordre 4 pas de temps 
adapt (ASRK4) et I'algorithme de Verlet. Ces algorithmes sont dcrits en dtail dans 
les paragraphes suivants. Les trois premiers ont un pas de temps fixe alors que le 
quatrime adapte le pas de temps en fonction de la variation du potentiel: si les vari- 
ations sont grandes, le pas de temps est petit; si les variations sont faibles, le pas de 
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2 Dynamique Moleculaire et Modele en Liaisons Fortes 



temps est grand. Un paramtre gre les variations du pas de temps, qui se traduisent 
par un cart dans I'espace des 6N coordonnes (3 d'espace et 3 de vitesse par particule). 

Dans une simulation de Dynamique Moleculaire, ctiaque atome est considere comme 
une masse ponctuelle, dont le mouvement est determine par I'ensemblc dcs forces ex- 
ercees sur lui par les autres atomes en fonction du temps. Le pas de temps doit etre 
approximativement 10 fois plus rapide que les mouvements les plus rapides que Ton 
veut simuler (Echelle de temps atomique: lfemtosec(10-15s). Grace aux puissants 
ordinateurs dont on dispose on pent simuler des trajectoires pouvant durer quelques 
nanosecondcs (10-9 secondes). Les forces exercees sent considerees comme constantes 
pendant le pas de temps. 

Les ordinateurs actuels peuvcnt traiter des systemes comportant jusqu'a des dizaines 
dc millions dc variables, cepcndant les tallies dcs systemes traites restent encore assez 
modestes en comparaison avec les ensembles macroscopiques. 

Ces nombres croitront probablement de maniere significative dans un proche avenir. 
Les mesures doivent aussi etre assez nombreuses pour atteindre une precision statis- 
tique suffisante. 



2.2.1 Methode d'Euler 

C'est la methode la plus simple. On cherche a resoudre un systeme de type: 

^ = fii{yjht) ie[l,n] (2.3) 

oil les fonctions yi{t) sont les inconnues et les fonctions fj({yj}, t) sont connues et 
dependent de la variable d'integration t mais aussi des valeurs prises par I'ensemble 

des Yj. 

II s'agit d'equations qui sont a priori du premier ordrc, alors que (1.1) est du deuxieme 
ordre. Ce n'est pas une limitation serieuse, car, imaginons que Ton ait une equation 
du type: 

il suffit d'introduire une variable supplementaire en posant: 

dz(t) 

"^'^ = -dT 

et Ton obtient un systeme de deux equations differentielles du premier ordre: 
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2.2 Simulation de Dynamique Moleculaire 



Ceci est generalisable a un nombre quelconque d 'equations du deuxieme ordre pour 
obtenir deux fois plus d'equations du premier ordre ^. Par exemple, imaginons un 
systeme de deux pendules simples couples. Les equations du mouvement d'un tel 
systeme s'ecrivent: 



1 = -fsin9i + -^(9 



72 



On peut fairc la dorrespondance: 



yi 


<— > 






<— > 






<— > 




y4 


<— ^ 




et: 






/i 




y2 


/2 




-f^sinyi + 


/3 






/4 







c 



mill 



^{ys-yi) 



c 



fiyi-ys) 



La methodc d'EuIcr, a proprement parler, consiste a remplacer le systeme d'equations 
(1.3) par un developpement au premier ordre: 

yi{t + dt) = yi{t) + fi{yj{t),t)dt 
oil dt prend maintenant une valeur petite mais finie. 

Le probleme ici est que ii{{yj{t)},t), la derivee, est evaluee a I'instant t, au debut 
de I'intervalle [t, t+dt], ce qui est un choix pour le moins arbitraire. La methodc de 
Runge-Kutta developpee ci-dessou tente de repondre a cette question en faisant une 
moyenne sur plusieurs evaluations. 

2.2.2 Methode de Runge-Kutta 

Developpement de la methode de Runge-Kutta a I'ordre 4, s'ecrit: on commence par 
calculer les derivees a I'instant t^ comme dans la methode d'Euler: 

i^ = M{yjm,t) 

puis on avance d'un demi-pas, et Ton reevalue les derivees: 



^et meme a un nombre quelconque n d'equations de degre p pour obtenir n x p equations du premier 
ordre 

^on s'interesse ici a des equations qui dependent du temps, mais evidemment, ceci est generalisable 
a n'importe quel type de variable 
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2 Dynamique Moleculaire et Modele en Liaisons Fortes 



On revient au point de depart, puis on fait de nouveau un demi-pas avec les nou- 
velles valeurs des derivees /]^', et enfin on evalue de nouveau les derivees: 

=..W + /f|/f = /.({yf(0},t + |) 

On revient de nouveau au point de depart et Ton fait un pas complet cette fois en 
utilisant la troisieme evaluation des derivees: 

L'on possede alors quatre evaluations des derivees dont on fait une moyenne ponderee: 

et Ton fait alors un pas complet avec ces "derivees": 

yi{t + dt) = yi{t) + fl^dt 

Un grand nombre de problemes de physique ont ete resolus de la sorte (d'autant 
plus que le programme correspondant tient en quelqucs ligncs), mais il faut se rap- 
peller qu'en general 99%, sinon plus, du temps de calcul dans une simulation de 
dynamique moleculaire est consacre au calcul des forces, c'est a dire les derivees en 
termes d'equations differentielles, a cause des "^"^"^^ paires d'atomes dont il faut cal- 
culer I'interaction. Ainsi, une methode dans laquelle il faut calculer les forces quatre 
fois par pas n'est probablcment pas la mcilleure, meme si l'on pent en attendre la 
possibilite de clioisir un pas d'integration plus grand. 

2.2.3 Algorithme de Verlet 

C'est un algorithme couramment utilise, tres simple celebre, appele encore "saute- 
mouton" ou " leap-frog " . L'idee de base est d'ecrire la fonction de position f{t) en 
troisieme ordre du developpement de Taylor 

Soit v{t) la Vitesse, et a{t) I'acceleration, et b la troisieme derivee de f{t) par 
rapport a t. 

En entree, on donne une configuration: 
{fi(0),....,rjv(0),i;i(0), ,%(0)} 

et on souhaite obtenir en sortie, apres un intervalle t et avec la meilleur precision 
possible, une configuration finale: 

{fi{t),....,rN{t),Vi{t), ,VNit)} 



22 



2.3 Necessite des methodes semi-empiriques: 



L'algorithme de Verlet, que I'oii se propose de justifier dans ce qui suit, permet de 
realiser cette operation de faon satisfaisante, en utilisant I'operateur de Liouville: 
f{t + 6t) = f{t) + v{t)6t + (l/2)a(t)<5t2 + {l/6)b{t)5t^ + 0{6t^) 
f{t - 6t) = r{t) - v{t)6t + {l/2)a{t)6t^ - {l/6)b{t)6t^ + 0{dt^) 
en sommant les deux equations on trouve I'equation d'evolution approximative 
pour une particule: 

r{t + 6t) = 2r{t) - f{t - 6t) + a{t)6t^ + 0{6t^) (2.4) 
I'acceleration est donnee par: 

a{t) = -(1/m) V V{r{t)) (2.5) 

et les vitesses : 

vit) = ri^ + ^t)-rit-St) ^^^^^ 
2dt 

La determination de la vitesse permet d'obtenir la position de I'atome a I'aide de 
I'equation precedemment enoncee a I'instant (t+5t). 

La repetition de cette procedure a des intervalles de temps discrets en fonction de 
la vitesse aboutit a 1' identification de la trajectoire. 



2.3 Necessite des methodes semi-empiriques: 

Les lois fondamentales necessaires a la comprehension des phenomenes en physique et 
chimie sont en grande partie connues. La difficulte est que leur application rigoureuse 
conduit souvent a des equations trop comDliauees|12j. II est done souhaitable de 
developper des methode pratiques approchees, en particulier dans le domaine de la 
mecanique quantique dans le but d'expliquer les particularites des systemes atom- 
iques complexes sans recourir a des calculs excessivements difficiles, de type "ab- 
initio" lesquels ne sont possible qu'avec un nombre limite d'atomes et necessite des 
calculateurs de plus en plus performants. Par exemple, la resolution des equations 
de Roothaan concernant les integrales de repulsion inter electronique de coulomb 
conduit a des calculs extremement volumineux. Par contre, les methodes approchees 
semi-empiriques simplifient les modeles en negligeant la plupart (ou la totalite) des 
integrales moleculaire de repulsion coulombienne. En outre, les integrales de coeur 
ne sont pas normalement calculees en toute rigueur, mais on les considere comme des 
parametres ajustables de fagon a obtenir la meilleur concordance avec des resultats 
experimentaux ou des calculs " ab-initio" . 
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2 Dynamique Moleculaire et Modele en Liaisons Fortes 



2.4 Modele en Liaisons Fortes: 

Dans le chapitre 2 nous allons utilise une mcthodc semi-empirique utilise le modele en 
Liaisons Fortes ("Tight-Binding"), pour I'etude du systeme CoPt. Cette approche, 
elle part d'une base que I'interaction electronique entre les atomes est supposee etre 
relativement faible et les electrons suffisament localises autour des noyaux, ce qui per- 
met d'ecrire leur fonction d'onde comme une combinaison lineaire d'orbitales atom- 
iques centrees sur chaque site. Cette approximation est souhaitable quand les dis- 
tances interatomiques sont grandes devant I'etenduc dcs fonctions d'ondes atomiques. 
Elle s'applique en particulicr au cas des bandes etroitcs (3d, 4d, 5d et 5f) des metaux 
de transition et a tout materiau caracterise par des etats de valence localises autour 
des atomes. 

Dans ce sous-chapitrc nous dccrivons I'expression de I'hamiltonien dans I'approximation 
des Liaisons Fortes. Puis nous en deduisons I'expression de la densite d'etats electronique 
locale en volume, qui permet le calcul de I'energie de bande par site. Puis nous ex- 
primons la densite electronique dans le cadre de I'approximation au second moment 
(SMA, "Second Moment Approximation ") qui permet d'obtenir un potentiel inter- 
atomique a " N-corps" . 

2.4.1 Hamiltonien en Liaisons Fortes: 

L'Hamiltonien a un electron s'ecrit: 



oil T est le terme d'energie cinetique et Vj le potentiel atomique centre sur I'atome 
i, chaque fonction d'onde pent s'exprinicr comme une combinaison lineaire d'orbitales 
atomiques notees \iX > oii i est I'indice du site atomique et A I'indice de degenerescence 
de I'orbitale {X = 1, . . . ,1). 

On se place dans une base supposee orthonormee et complete des orbitales atomiques 
(ce qui necessite que le recouvrement des orbitales soit negligeable) de telle sorte que: 



N 




(2.7) 



i=l 



< iX\jn >= Sijdx, 



(2.8) 



Dans cette base, les elements de matrice de I'hamiltonien peuvent s'ecrire: 




< iX\T + ^ Vfc|j> >=< iX\T + Vj\jfi > + <iX\J2 Vk\jl^ > (2.9) 
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2.4 Modele en Liaisons Fortes: 



Sachant que chaque orbitale atomique obeit necessairement a I'equation de Schrodinger 
pour un atome isole et dans la mesure ou Yj est assimilable au potentiel d'un atome 
libre, on a: 

(r + V,)|j>>=ej^|i/i> (2.10) 

ou ejfj, est I'energie du niveau atomique de I'atome j pour I'orbitale qui peut 
dependre de la structure atomique locale. Ceci permet d'ecrire: 

H-f = ejf,5,j6x^ + ^ < iX\Vk\jfi > (2.11) 

On peut remarquer que le second terme est constitue d'integrales a trois centres et 
d'integrales a deux centres. Les premieres sont generalement negligees et les secondes 
peuvent se limiter aux proches voisins du fait que les orbitales d s'attenuent exponen- 
tiellement avec la distance interatomique. Ainsi, parmi les elements < iX\Vk\jiJ- >, 
on ne conserve que ceux ou k=i. Les elements de matrice intra-atomiques sont 
representes par: 

H-l" = ejM'^AM+ < J2 ^k\m > (2.12) 

ovL le second terme, connu sous le nom d'integrale de champ cristallin ou integrale 
de derive, a pour effet d'abaisser legerement le niveau atomique. II est habituellement 
neglige. 

Les elements de matrice inter-atomiques appeles integrates de saut ou de transfert 
sont de la forme: 

Hlf=<tXmj,j.>=pl^ (2.13) 

Elles sont responsables de la formation de la bande a partir des niveaux atomiques 
discrets (i.e. levee de degenerescence atomique). En d'autre termes, elles permettent 
aux electrons de "sauter" de site en site dans le solide. Elles s'expriment en fonc- 
tion des cosinus directeurs 1, m, n du vecteur Rij=Ri+Rj et d'un certain nombre 
d'integrales independantes appelees parametres de Slater et Koster pTSj . 

Cela nous permet de formuler I'hamiltonien en Liaisons Fortes en fonction des 
parametres Cix (energie du niveau atomique) et P^j^ (integrale de saut) de la maniere 
suivante: 

H = ^\iX>ea<iX\+ \i^>Ptf<jf^\ (2-14) 

i,X i,j^i,\,\ytfi 

OU les sites atomiques i et j sont proches voisins. 
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2 Dynamique Moleculaire et Modele en Liaisons Fortes 



Selon la methode usuelle, on resout I'equation de Schrodinger: 



H\n >= En\n > (2.15) 

ou les etats propres \n > sont des combinaisons lineaires d'orbitales atomiques, 
soit: 



\n >=^^Ci\\iX > avec I ^_^]^ ' "'^ (2-16) 



Ce qui conduit a un systeme lineaire de dimension NxZxNx/ que Ton reduit a 
Ix I si on considers un solide avec 1 atome par maille et dans lequel on applique le 
theoreme de Bloch: 



cu = -^e^'-'^'c^ik) (2.17) 



Les equation du systeme Ix I sont alors de la forme: 



{E^ - E)c^{k) + J2 (ixi^e'^-'^'c^ik) = (2.18) 

Leurs solutions en tout point k de la zone de Brillouin caracterisent la courbe de 
dispersion E(k) et la densite d'etats resultante. Ceci permet egalement d'obtenir 
I'energie de bande. Cependant, cette methode necessite, pour utiliser le theoreme de 
Bloch, d'avoir des conditions periodiques dans toutes les directions. 



L'avantage du formalisme des Liaisons Fortes est justement de pouvoir caracteriser 
la densite d'etats dans I'espace direct, sans avoir a diagonaliser I'hamiltonien. 



2.4.1.1 Densite d'etats electroniques: 

Tres souvent, on suppose que I'alliage pent etre decrit par une structure de bandes 
dont la densite d'etats n(e, c) correspondant a une concentration c se deduit simple- 
ment de la densite d'etats du metat pur{n(£, 0) = n{e)}, sans changement de forme, 
par un simple deplacement des energies: 

La densite d'etats electronique totale n(E) est simplement le nombre d'etats d'energie 
compris entre E et E + dE, divise par dE. On pent I'ecrire en utilisant la notation 
de Dirac: 

<E) = ^Y.^{E - E^) (2.19) 

n 
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2.4 Modele en Liaisons Fortes: 



Si on definit I'operateur 5{E.Id — H), ou Id est la matrice identite, tel que: 

5{E.Id-H)\n>=5{E-En)\n> (2.20) 
alors on peut exprimer la densite d'etats comme: 

n{E) = -^Tr6{E.Id - H) (2.21) 

ou la trace est definie sur I'ensemble des etats electroniques de I'espace de Hilbert 
pour une base quelconque, par exemple la base des orbitales atomiques \iX > decrite 
precedemment . 

On peut definir ainsi la densite d'etats locale nj(E) au site i a partir des densites 
d'etats projetees sur chaque orbitale A notees ni;^(E) [Hj : 

ni{E) = jJ2x = llnix{E) (2.22) 
sachant que les nj;s^(E) sont definies comme: 

na{E) = ^ nc*^{En)ca{E„)6{E - E^) (2.23) 
telles que d'apres la condition de normalisation: 

/+00 r+oo 
nix{E)dE = / ni{E)dE = 1 (2.24) 
-oo J ~oo 

La densite d'etats totale est alors definie par la somme normalisee de toutes les 
densites locales: 



N IN 

= nY.''^^^) = JnT^Y. "'^(i?) (2.25) 

i=l A=l i=l 

II existe deux moyens de caracteriser la densite d'etats electroniques 6j la methode 
des moments qui permet de definir n(E) par la donnee de ses moments d'ordre p, 
p etant plus ou moins eleve selon le degre de precision souhaite, et la methode de 
recursion utilisant la fraction continue [Hj- Ces methodes sont reliees par des relations 
simples (du moins en ce qui concerne les premiers moments). 

2.4.1.2 Methode des moments 

D'une maniere generale, la densite d'etats peut etre caracterisee par la donnee de 
tons ses moments, le moment d'ordre p s'ecrivant: 

r+co 

Hp= EPn{E)dE. (2.26) 
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D'apres (fT^ . et en utilisant le fait que E'^Tr5{E.Id - H) = TrH'^5{E.Id - H) 
jl4j . on en deduit que le moment d'ordre p s'exprime egalement par: 

^Ji, = ^TrH^ = E < ^^l^'l^^ > (2.27) 

iX 

La deuxieme egalite est obtenue en developpant la trace dans la base des orbitales 
atomiques. Cette notation est particulierement interessante puisqu'en introduisant p 
relations de fermeture X^j/x liA* iA*l=l on obtient p equations de la forme: 

/ip = ^ E <iMH\jfJ-><jfJ-\H\ka>---<...\H\iX> (2.28) 

iX,j ii,ka,... 

D'apres I'expression de H donnee par (|2.14|) . on obtient une relation simple entre 
le moment d'ordre p et les parametres e?^ et P^j^. Cette relation revient a considerer 
tons les chemins fermes de p sauts partant du site i (et de I'orbitale A) et y revenant, 
y compris les sauts "sur place" concernant les niveaux atomiques. 

Les premiers moments se calculent aisement et out une signification physique: 
Le moment d'ordre zero represente le nombre total d'etats disponibles (norme) 

iX 

le moment d'ordre un correspond au centre de gravite de la bande 

^1 = 7^7 E < ^^I^K^ >= 7^ E ^°A, (2.30) 



IN ^ ' ' IN 

iX iX 



le second moment: 



iX i,jj^i,X,iiytX 

dont la racine carree represente I'ecart quadratique moyen de la densite detats 
(dans le cas d'une bande d seule, la racine carree du second moment pent etre reliee 
a la largeur a mi-hauteur de la bande d). 

/i3 renseigne sur I'asymetrie de la bande, //4 sur son etalement. . . . 

La donnee de tons les moments permet de d'ecrire la densite d'etats grace a la 
fonction caracteristique des moments: 



n=0 



/W = E^^/^« (2.32) 
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2.4 Modele en Liaisons Fortes: 



dont la transformee de fourier donne la densite d'etats jl5j : 

1 f + oo 

<E) = ^j^ e'^^'frn. (2.33) 

Comme 11 est impossible de connaitre tous les moments, on pent approximer la 
densite d'etats par une fonction ayant les meme premiers moments, ou bien on 
utilise la fraction continue qui permet de calculer directement la densite d'etats par 
I'intermediaire de coefficients que Ton peut relier aux moments. 



2.4.1.3 Expression de I'energie de bande 

L'energie de bande au site i s'exprime par: 

= E r^^^ - ^ixHxiE,e,x)dE (2.34) 
avec, d'apres ce qui precede pour les sites de surface: 

= 6° + 6e^x (2.35) 

oil e° est le niveau atomique de I'orbitale A dans le volume et 5eix est le decalage 
de ce niveau correspondant au site i. D'apres ce qui precede, 6€i\ 7^ lorsque i est 
un site de surface. 

Si on developpe I'expression (|2.34|) . on obtient: 

^ f-Ep 

Ei=J2 - ^x)^i^(E, e,x)dE - NixSax (2.36) 

A -^-"^ 

Oil Nix est le remplissage de la bande A au site i. 

Dans le cas d'une neutralite globale (deplacements rigides des niveaux atomiques 
de toutes les orbitales) on a: 

5eix = 6ei VA et N^a / N^x mais Ea ^ix = 

ou est le nombre total d'electrons sur toutes les orbitales et done: 

Y,^ixS^iX = N''6ei (2.37) 

X 

Dans le cas d'une neutralite par orbitale (pas de transfer de charge entre les or- 
bitales): (5e.jA / Sei^ si A / et Nj^ = V i et done: 

Y,NixSeix = N°x^eix (2.38) 

A 
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Dans les expressions ci-dessus, la densite d'etats local nj;^(E, 5ei\) est celle determinee 
apres le calcul autocoherent tel qu'on I'a explicite precedemment. On pent egalement 
determiner I'energie de bande de maniere raisonnable a partir de la densite d'etats 
initiale (avant autocoherence) notee n*^(E) en autorisant un deplacement du niveau 
de Fermi pour conserver la charge de volume. L'energie de bande est alors ap- 
proximee par: 

^'^E r {E-el)nUE)dE 

^ J-mfty 

En introduisant le niveau de Fermi du volume selon la formule 

^ E( - '>txiE)dE - {E - el)nUE)dE) (2.40) 

on peut approximer (|2.39|) dans le cas oii E^ est voisin de E^? par I'expression 
suivante: 

^ E( Ev*.xiE)dE - EF5Ni^ - e^N^^) (2.41) 

Pour determiner I'energie de site totale, il faut ajouter au terme de bande attractif 
ci-dessus un terme repulsif correspondant a la repulsion entre les noyaux atomiques. 
Celui-ci est determine empiriquement par un potentiel de type Born-Mayer: 

= ^e-P(^~^) + C (2.42) 
j 

oil A, p et C sont ajustes sur les valeurs experimentales de I'energie de cohesion, 
du parametre de maille jTj et du module de compressibilite [T^. II faut noter que 
la constante C permet de compenser I'erreur faite sur les niveaux atomiques in- 
tervenant dans le terme attractif. En effet, les niveaux atomiques etant issus d'un 
schema d'interpolation ^3] ne sont qu'effectifs et ne peuvent en aucun cas rendre 
compte des niveaux atomiques de I'atome isole. 

2.5 Approximation au Second Moment de la densite d'etats 
(SMA) 

La cohesion des metaux de transition est essentiellement gouvernee par les electrons 
d ^j. De plus, J. Friedel^Tj et F. Ducastellej^O ont montre qu'il n'est pas necessaire 
de connaitre les details de la densite d'etats pour caracteriser les proprietes de 
cohesion de ces metaux et leurs evolutions le long de la serie lorsqu'on fait varier 
le nombre d'electrons d. II suffit d'avoir la repartition globale en energie de tons les 



(2.39) 
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2.5 Approximation au Second Moment de la densite d'etats (SMA) 



n(E) 



lOAV 



E 



-W/2 



+W/2 



Figure 2.1: Densite d'etats approximee au second moment. 



etats, qui est fournie precisement par la largeur moyenne de la bande. L'approximation 
au second moment consiste a caracteriser la densite d'etats en se limitant aux mo- 
ments d'ordre inferieur ou egal a 2 que nous avons developpes precedemment dans 
le cas general (|2.29l a \2.'A1^ et que nous utilisons a present en ce qui concerne les 
orbitales d (1=10) pour le site i. 

Si on choisit une forme rectangulaire pour nj(E), centree en ef choisie comme 
origine des energies et telle que figure (|2.1|l : 



(on remarque qu'on a choisi ici de normer la densite a 10, nombre d'etats disponibles 
en bande d). L' identification du second moment de cette densite d'etats de forme 
rectangulaire est triviale et conduit a: 



On en deduit que nj(E)=10-v/l2^ • /?• . Sachant que le nombre d'electrons par 



G [-Wi/2,+Wi/2],Wini{E) = 10 



(2.43) 




(2.44) 
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atome dans le metal est donne par: 

/Ep 
ni{E)dE, (2.45) 
-oo 

on pent calculer le niveau de Fermi: 

Ef = {^- \)W. (2.46) 
et ainsi I'energie de bande au site i d'apres I'expression: 

E\ = 12 En,{E)dE = -N,{1 - ^)^. (2.47) 
D'apres 1)2. 44() et si on pose de plus: 

A^/ = /?o^e''^^~'^ete = l3l^V?>Nd{lQ - N^) (2.48) 

ou Tij est la distance entre les sites i et j et ro la distance d'equilibre dans le volume 
entre atomes premiers voisins, alors I'energie de bande au site i s'ecrit de la fagon 
suivante: 

Et = -^ery^^' (2-49) 

Notons que I'application de ce modele aux metaux nobles peut sembler paradoxale 
puisque ces metaux ont une bande d pleine, N(i=10 et done ^=E^=0! Pour expliquer le 
succes constate empiriquement de la formule (|2.49() Dour ces metaux, B. Legrand et M. 
Guillope ^Hl ont montre que dans ce cas, la cohesion n'est plus due a I'elargissement 
de la bande, mais a son deplacement en fonction de I'environnement local. On montre 
effectivement que I'integrale de derive correspondant au deplacement de la bande, 
n'est plus negligeable devant I'elargissement de celle-ci (I'integrale de saut) pour les 
metaux de fin de serie de transition, done a fortiori pour les metaux nobles. Cela est 
probablement lie a I'importance de I'hybridation s-d en fin de serie. 

L'energie de bande exprimee par ()2.49|) a partir de la structure electrique va con- 
stituer le terme attractif du potentiel en Liaisons Fortes dans I'approximation au 
second moment. On peut remarquer, a ce stade, que les parametres du modele (^ et 
q) sont ajustes sur la structure electronique. 

Finalement, I'energie par atome est obtenue en ajoutant au terme attractif de bande 
1)2. 49(1 un terme repulsif sous forme d'interactions de paires, pour assurer la stabilite 
du reseau, soit: 

Ei = E\ + El (2.50) 
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L'energie au site i s'exprime alors sous la forme: 



^^ = -J E ejjexpi-2qui^-l))+ Ajjexp{{-pjj{^-l))) (2.51) 

y j,rij<rc ^ j,rij<rc ^ 

oil (^,q,A et p) sont les parametres du modele (la constante C est nulle dans ce cas 
ou Ton traite uniquement les electrons d). D'apres ce qui precede, ^ et q pourraient 
etre determines directement a partir de la structure de bande dans le cas oii on ne 
tient compte que de la bande d. 

Les prametres du potentiel son ainsi ajustes pour reproduire les valeurs experiment ales 
de l'energie de cohesion, du parametre de reseaujlj et de quelque constantes elastiques 
[H] (module de compressibilite B et modules de cisaillement C44 et C ) du metal. II 
est essentiel de remarquer que la seule presence de la racine carre dans (|2.49j) confere 
au potentiel un caractetre non additif, autrement dit " a N - corps" , primordial dans 
I'obtention de quelques resultats remarquables que les potentiels de paires sont inca- 
pables de reproduire, par exemple: 

- Les energies de formation de lacune jl9j. 

- Les relaxations vers I'interieur (contraction) des surfaces \2{)\ I18| 
Actuellement, seule les methodes ab initio (LDA, LMTO, car-Parnello) issues de 

la theorie de la fonctionnelle de la densite calculent l'energie totale du systeme sans 
parametrisation. Cependant elles sont relativement lourdes a mettre en oeuvre pour 
determiner la structure et la morphologie dequilibre d'agregats de metaux de transi- 
tion, meme au niveau des plus petites failles. 

Une approche semi-empirique comme le potentiel en Liaisons Fortes dans I'approximation 
au second moment constitue une methode de choix a I'heure actuelle pour caracteriser 
de maniere realiste des particules de faille comparable a celle des particules observees 
experimentalement, c'est a dire de plusieurs centaines a plusieurs milliers d'afomes. 

2.6 Methode Monte- Carlo 

Les simulations Monte Carlo consistent en une marche dirigee dans I'espace des con- 
figurations possibles du systeme. L'algorithme est congu de maniere ce que chaque 
configuration apparaisse au cours de la simulation avec une probabilite egale a son 
poids statistique. Les moyennes thermodynamiques sont par consequent egales, dans 
la limite d'un nombre infini d'iterations Monte Carlo, aux moyennes realisees sur 
les configurations generees lors de la simulation. Elle porte ce nom parce qu'elle est 
basee sur I'utilisation de nombres aleatoires. 

Les travaux les plus anciens font etat des premiers succes qui ont etabli la simulation 
comme un outil majeur d' investigation de la physique de la matiere condensee. La 
methode de Monte-Carlo (MC) fut developpee par Von Neuman, Ulam et Metropo- 
lis, a la fin de la seconde guerre mondiale, pour I'etude de la diffusion des neutrons 
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dans un materiau fissile. N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. 
H. Teller et E. Teller furent les pionniers de I'investigation de la matiere par 
simulation sur I'ordinateur. lis realiserent la simulation d'un liquide simple (disques 
durs se deplagant en deux dimensions) par la methode MC. lis proposerent ce qui 
porte desormais le nom de MC Metropolis et qui est devenu la base des simulations 
MC des systemes de particules en interaction. 

Cette methode permet I'estimation des moyennes de grandeurs physiques donnees 
par la formulation de Gibbs de la mecanique statistique sous la forme d'integrales mul- 
tidimensionnelles. Les premieres simulations furent realisees dans I'ensemble canon- 
ique (N, V et T constants), puis la technique fut etendue aux autres ensembles statis- 
tiques. On genere une sequence aleatoire d'etats accessibles (chaine de Markov) dans 
I'espace des configurations du systeme. On echantillonne en privilegiant les regions 
ou le facteur de Boltzmann (exp(-U/KBT)), c'est-a-dire la densite de probabilite 
de I'ensemble canonique dans cet espace, est le plus eleve (algorithme de Metropo- 
lis). La probabilite d'une configuration particuliere d'energie potentielle Uj est alors 
proportionnelle a exp(-Ui/KBT). 

2.7 Comparaison entre methode MC et DM 

Les potentiels inter-atomiques utilises dans la dynamique moleculaire classique souf- 
frent souvent de la faon tres approchee et indirecte dont les electrons sont traites. 
A priori, I'effet des electrons est cache dans les parametres du potentiel, qui doivent 
etre determines par ajustement. En general, cet ajustement se fait sur des donnees 
experiment ales, qui ne sont pas toujours facilement accessibles, ou par des resultats 
de calculs " ab-initio" . 

La choix entre la methode Monte-Carlo ou la methode Dynamique Moleculaire 
se fait en considerant I'efficacite relative de ces deux techniques en fonction du 
probleme a traiter. Pour I'obtention des proprietes statiques d'equilibre des systemes 
(quasi-)ergodiques, les deux methode sont equivalentes, car dans ce cas les moyennes 
temporelles et les moyennes d'ensemble fournissent les mmes resultats. II est alors 
preferable de choisir I'approche DM puisqu'elle fournit le mouvement des particules. 
En effet, elle permet d'atteindre en outre les proprietes dynamiques et de trans- 
port qui sont inaccessibles par la technique MC. Le mouvement detaille des atomes 
pent tre analyse et conduire notamment a la connaissance des mecanismes precis des 
phenomenes physiques (diffusions dans les solides). 

D'autre part, la DM offre un echantillonnage plus efficace de I'espace des con- 
figurations quand il s'agit de traiter des transitions de phases structurales ou des 
changements conformationnels de grandes molecules. Dans ce cas, le chemin a par- 
courir d'une region de I'espace des phases a une autre requiert un rearrangement 
collectif des coordonnees de nombreuses particules. La DM permet alors de trouver 
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2. 7 Comparaison entre methode MC et DM 



des chemins plus directs que ceux generes par des mouvements MC aleatoires non 
correles et de ce fait moins efficaces. 

Par centre, la methode MC est plus facile a appliquer, et surtout plus interessante, 
dans le cas d'une energie potentielle mathematiquement compliquee, puisqu'on pent 
s'affranchir d'une evaluation explicite des forces, difficile mais indispensable a la DM. 
Par ailleurs, une situation oii la nature non rccUc des mouvements MC pent etre 
exploitee est celle des melanges, notamment quand I'interdiffusion des especes est 
trop lente pour tre observee a I'echelle de temps de la DM. Enfin, de nombreuses 
techniques MC ont ete developpees specialement pour le calcul des energies libres. 
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2 Dynamique Moleculaire et Modele en Liaisons Fortes 



36 



3 Modele TB-SMA pour I'etude du CoPt 



3.1 Introduction 



Ce chapitre commence d'abord par un rappel des principales proprietes de ce systeme, 
nous decrirons ensuite le modele semi-empirique TB-SMA et les methodes permet- 
tant la parametrisation de ce modele. Cette etude presentee a ete entreprise dans le 
but d'estime la validit d'un potentiel second moment dans I'extension de son utilisa- 
tion aux alliages CoPt, la methode que nous avons utilise consiste a reproduire des 
donnees experimentales de I'alliage equiatomique CoPt (LIq). Par contre une etude 
publiee[22] a utilise une approche differente en utilisant des donnees experimentales 
de I'alliage CoPt tres dilue, pour des tudes de couches minces dposes sur un substrat 
pur Pt ou pur Co. Dans tout le memoire, nous utiliserons les potentiels obtenus 
par les deux methodes notees Potentiel 1 et Potentiel 2, nous comparons les deux 
approches en referent. 

Nous exposons finalement les resultats obtenus par Dynamique Moleculaire concernant 
les barrieres energetiques. 



37 



3 Modele TB-SMA pour I'etude du CoPt 



3.2 Proprietes du systeme Co-Pt 

Parmi tous les alliages deja etudies, le systeme CoPt est Tun des plus prometteurs 
car susceptible de developper de fortes energies d'anisotropie magnetocristalline du 
fait de la stabilite de la phase anisotrope LIq et du couplage spin-orbite tres eleve 
dans le compose equiatomique ordonne de structure LIq. Ceci a ete mis en evidence 
par Brissonneau et al. , et confirme par Eurin |H l^j . 

Ce systeme a fait I'objet de tres nombreux travaux tant experimentaux que theoriques, 
et ce sont ses proprietes catalytiques et magnetiques qui font son interet. Nous 
presenterons ici ses proprietes structurales et cinetiques (migration atomique), ainsi 
que ses proprietes magnetiques. 

3.2.1 Diagramme de phases 

Le diagramme de phases du systeme CoPt donne par Hansen est relativement 
imprecis: il mentionne I'existence des phases ordonnees LIq et LI2 autour des com- 
positions CoPt et CoPts respectivement. Les contours du diagramme de phases 
chimique et magnetique du systeme Co-Pt ont ete precises par Dahmani |24| I25j en 
utilisant la diffraction de rayons X sur des echantillons trempes, et des mesures de 
resistivite in situ. Ce diagramme a ete complete par C. Leroux [22j en utilisant la 
microscopic electronique a transmission (MET) in situ sur des echantillons trempes. 

La figure H3.2.1|) presente un apergu general de ce systeme, ainsi que les temperature 
de curie des phases ordonnees et desordonnees. Notre etude concerne exclusivement 
la composition equiatomique 50-50 qui fait partie du domaine d'existence de la phase 
LIq. Cette phase ordonnee est constituee d'une succession de plans purs cobalt et 
purs platine dans la direction [001] (figure 2.2). 

Comme une grande partie de ce travail est consacree a la croissance et aux pro- 
prietes de I'alliage CoPt de structure LIq, il est important de mentioner quelques pro- 
prietes structurales de cette phase. La transition ordre-desordre qui se situe a 1110 K 
a temperature croissante figure (2.3 a) et a 1070 K a temperature decroissante [22 est 
fortement du premier ordre, c'est-a-dire que la discontinuite du parametre d'ordre a 
la temperature de transition Tc est grande. La mise en ordre a longue distance dans 
cette structure s'accompagne d'un changement de symetrie cristalline d'une structure 
tetragonale (tfc) a cubique a face centrees (cfc). Par rapport a la maille cfc de la 
phase desordonnee, la maille tfc correspond a une reduction du parametre de reseau 
suivant I'axe [001] (axe "c") et a une dilatation suivant les directions [100] et [010] 
(axes "a") figure(2.3 b). Des simulations du diagramme de phase par la methode de 
variation des amas f CVM) [?7| ont montre que le parametre d'ordre varie de 1 a 0.85 
entre et Tc figure (2.3 a), ce qui signifie que I'ordre dans les domaines ordonnes est 
eleve, des I'apparition de la phase ordonnee, et la discontinuite du parametre d'ordre 
a la temperature de transition Tc est grande (il varie brutalement de 0.85 a 0, voir 
figure (|3.2.1l al). Par consequent, dans un alliage massif a I'equilibre, un parametre 
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3.2 Proprietes du systeme Co-Pt 




Figure 3.1: Diagrame de phase CoPt 



d'ordre inferieur a 0.85 indique que I'alliage est constitue d'un melange de phase or- 
donnee et desordonnee. Tous ces resultats theoriques ont ete confirme indirectement 
par I'etude in situ des parametres de reseau par diffraction des rayons X en fonction 
de la temperature |2S| figure(2.3 b) qui a montre qu'une forte correlation entre le 
parametre d'ordre et les parametres de reseau. 

Au sein de notre laboratoire, le systeme FePd a ete etudie par resistivite 'in-situ' et 
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3 Modele TB-SMA pour I'etude du CoPt 
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Figure 3.2: Structure LIq autour de la composition AB et environement des deux 
atomes. 



RX par Messad et al. 
I'ordre de I'electron-volt. 



Les energies de migration trouvees sont plus faibles, de 



3.2.2 Mobilite atomique 

L'energie d'activation pour la mise en ordre (relaxation ordre-ordre) qui represente 
en premiere approximation, la somme des energies de formation et de migration des 
atomes a ete determinee par resistivite in situ |24j a basse temperature dans la phase 
CoPts- La valeur trouvee, 3.12 eV est en accord avec celle deduite de la croissance 
des domaines ordonnes etudiee par diffraction X [SHI et elle est proche des energies 
d'activation d'auto-diffusion dans le Cobalt et le Platine. 



3.2.3 Proprietes magnetiques 

Dans I'alliage CoPt les deux elements sont magnetiques, et I'alliage est un materiau 
ferromagnetique. En ce qui concerne le lien entre I'etat d'ordre et les proprietes 
magnetiques, des etudes ont montre que dans les alliages massifs CoPt de 

composition equiatomique I'aimantation a saturation est peu sensible a I'etat d'ordre 
(M°^'^/M^^*°'''^=1.06), tandis que la temperature de Curie diminue notablement avec 
le parametre d'ordre (Tc°'''^/Tc'^^^'""'^=0.88). 
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3.2 Proprietes du systeme Co-Pt 




I. ^ I 
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Figure 3.3: Proprietes du systeme CoPt, (a): Evolution du parametre d'ordre avec 
la temperature dans I'alliage CoPt massif, calculee par CVMJ'27', (b): 
Evolution des parametres de maille avec la temperature da I'alliage CoBl 
ordonne (a*^, c) et desordonne (a^)[22], (c): Effet de la composition sur 
la tetragonalite de la maille parametre de maille dans le plan (□) et 
perpendiculaire au plan (A), et rapport a/c (0). 



3 Modele TB-SMA pour I'etude du CoPt 



La repartition des moments magnetiques entre le cobalt et le platine reste un 
sujet de discussion. D'apres les mesures de diffusion magnetique des neutrons dans 
CoPts et CoPt, effectuees respectivement sur monocristal jHO] et polycristal |^ . 
les moments du cobalt et du platine seraient sensiblement les mmes dans les deux 
composes, et il semblerait d'apres ces resultats, que dans le systeme CoPt les moments 
du Co et du Pt varient peu avec la concentration, indiquant une faible dependance 
des moments avec leur environnement. 

Une estimation de I'aimantation a saturation a ete faite par Eurin,4 a partir des 
mesures effectuees sur un compose ordonne CoPt a temperature ambiante: I'etude 
de Dahmani on montre qu'il faut appliquer un champ d'environ 15 T pour sat- 
urer un polycristal ordonne. Par consequent, I'energie d'anisotropie magnetique d'un 
monocristal qui comporte un seul variant (en proportion de 97%), determinee par 
Eurin a partir de mesures effectuees a I'ambiante avec un champ magnetique maxi- 
mum de 7 T, est une limite inferieure. Cependant, sa valeur (46x10^^ erg/cm'^) est 8 
fois superieur a celle du Co dans la structure hexagonale compacte. 



3.3 Dynamique Moleculaire en Liaisons Fortes 
(TB-MD-SMA) 

L'algorithme de Dynamique Moleculaire Trempee en Liaisons Fortes avec I'approximation 
au second moment de la densite d'etats (SMA), permet de decrire le mouvement (les 
trajectoires) au cours du temps des atomes i en interaction a partir de positions 
hors-equilibre r^, soumis a des forces de rappels Fi vers leurs position d'equilibre par 
integration des equations du mouvement: 

- dvi d^r, 

dt dt^ 

Fi est la force agissant a I'instant t sur I'atome i de masse m, 
Vi est la Vitesse de I'atome i au meme instant 

fi est la position i a I'instant t, qui est obtenue a chaque instant par 
I'algoriyhme de Verlet |S2] 

Les coordonnees (Qj=Xj, y^, Zj) d'un atome i a un instant s'ecrivent: 

pa 

aiit + dt) = 2xi(t) -aiit- dt) + -^dt^ + eidf^). 

m 

La procedure de trempe consiste a "refroidir" le systeme en annulant la composante 
de la Vitesse d'un atome i des que: vf < 0, afin de le forcer a converger 
vers la position d'equilibre a temperature nulle. Cette procedure conduit done a la 
minimisation de I'energie potentielle a nK|33j. La force s'exergant sur chaque atome 



(3.1) 
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3.3 Dynamique Moleculaire en Liaisons Fortes (TB-MD-SMA) 



i, due a ses interactions avec les autres atomes, derive de I'energie potentielle du 
systeme de N atomes: 

Fi = — avec Etot = i-iEi 

an 

ou I'energie par site Ej est definie comme la somme d'un terme attractif de bande 
et d'un terme repulsif, comme vu au chapitre 1. 
Pour rappel, I'energie d'un site s'ecrit: 



^^ = -J E euexp{-2qjji^-l))+ Ajjexpii-pui^-m (3.2) 

y j,rij<rc ^ j,rij<rc 

Oil 



I, J: precise la nature chimique des atomes sur les sites (i,j) I, J G {Co,Pt}. 
^/j (Integrale de saut effectives) et qij: parametres de terme de bande. 
A/j et p/j parametres du terme repulsif. 

T^j correspond a la distance d'equilibre entre premiers voisins dans le 
metal pur I. 

^IJ— 2 

Vij est la distance entre les sites i et j. 

Vc est le rayon de coupure des interactions que Ton prend egal a la distance 
des seconds voisins. Audela de cette distance, les interactions s'ecrivent 
sous forme d'un polynome du cinquieme degre s'annulant de fagon con- 
tinue et derivable a la distance des troisiemes voisins, ceci pour eviter les 
effets lies a une discontinuite de potentiel. 



3.3.1 Parametrage du potentiel SMA 

Le potentiel SMA s'ecrit en fonction de differents types de parametres correspondant 
aux interactions Co-Co, Pt-Pt et Co-Pt (interaction mixte) ce qui fait en tout 12 
parametres a determiner. 

Parametres Co-Co et Pt-Pt: les valeurs des parametres (C/ZiAz/jPz/jq//) sont 
ajustes pour reproduire I'energie de cohesion, le parametre de maille, le module de 
compressibilite et I'equation universelle |35 des metaux purs. 

L 'equation universelle est une relation obtenue empiriquement, meme si elle a regu 
par la suite quelques justifications dans le cadre de I'approximation des Liaisons 
Fortes [H]. 

Elle relie I'energie du systeme, judicieusement normalisee, au parametre de reseau, 
lui aussi correctement normalise. 



43 



3 Modele TB-SMA pour I'etude du CoPt 



La relation obtenue s'est revelee valable pour une grande classe de materiaux 
(molecules, metaux purs, alliages metalliques, oxydes,. . . ). Cette relation est la suiv- 
ante: 

E*{R*) = -{l + R* + 0,05R*')e-^* (3.3) 

oil 

E* et R* sont sans dimension et sont relies respectivement a I'energie totale, E(R), 
et a la distance interatomique, R, [HSl par: 

E*{R*) = ^ 

"R* R—Rq 

^ - \Ro 

Ou 

A est une longueur caracteristique sans dimension donnee par A=(g^)^/^ 
B est le module de compressibilite 
Lo est le volume atomique 
Ec est I'energie de cohesion 

Ro est la distance d'equilibre entre premiers voisins a pression nulle. 

Nous definissons egalement les derivees premieres et secondes de I'expression (2.2): 

^= (R*-0,15 R*' + 0,05 R*' )e-^* 
^= (1-1,3R* +0,3 R*' + 0,05 R*' )e-^* 



Ainsi, lorsque la distance R est la distance d'equilibre Rq (i.e. a pression nulle), on 
obtien: 

E*(i?*) = -1 
/ dE* \ _ n 

[-^^)r=Ro - i 



Get ajustement permet d'obtenir les parametres du modele SMA pour les metaux 
purs Pt et Co, leurs valeurs sont donne dans le tableau (2.1). 

Parametres Co-Pt: Les parametres croise (^/j,A/j,p/j,q/j pour I / J) car- 
acterisent les interactions mixtes Pt-Co. II y a plusieurs methodes de les ajuster. 

i) Goyhenex et al.^\ ont ajuste les paraetres Go-Pt sur I'energie de dissolution 
d'une impurete Go dans une matrice (Pt) et de Pt dans (Go) et ont verifie qu'on 
se rapproche du paramtre de reseau exprimental, en calculant I'nergie de la configu- 
ration LIO en fonction de la distance, le potentiel ainsi obtenu est appele potentiel 
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3.4 Energie d'une configuration chimique {pf } 



Par. 


Pt-Pt 


Co-Co 


An 


0.242 eV 


0,189 eV 




2,506 eV 


1,907 eV 


VII 


11.14 


8.8 


q// 


3.68 


2.96 



Table 3.1: Parameters du potentiel SMA pour les metaux purs Co et Pt ajustes 
sur les valeurs experimentales (a, Ecoh et B) et a I'aide de I'equation 
universelle[2^. 

2. L 'inconvenient de cette demarche est que les energies de dissolution peuvent etre 
considerees comme une propriete de I'alliage CoajPti-^^ tres dilue, alors que notre but 
est de modeliser I'alliage equiatomique LIq. 

ii) Notre approche consiste a ajuster les parametres Co-Pt afin de reproduire 
raisonablement les proprietes (energie de formation et distance interatomique) de 
I'alliage equiatomique CoPt LIq. Notre potentiel est appele Potentiel 1, et son 
parametrage est expose en 2.5. 

Nous definissons ci dessous les energies de formation et de dissolution necessaire 
pour le parametrage du modele. 

3.4 Energie d'une configuration chimique {p"} 

L'expression de I'energie de la configuration chimique notee {pf} qui represente 
I'ensemble des facteurs d 'occupation de site pf tels que pf =1 si le site i est oc- 
cupe par un atome de type a (q=A,B) et pf =0 sinon. En utilisant l'expression du 
potentiel (2.2) I'energie de cette configuration s'ecrit alors: 

^(K}) = E E pf(-J H E ^'i^/V^p(-29/^(^ - 1))+ 

i a=A,B y j,rij<rcb=A,B '^O 

E E pj^ij^Mi-Pui^ - m (3.4) 

oil Iq'^ est la distance interatomique de reference utilisee pour le calcul des integrales 
de saut premiers voisins (nous choisirons ici celle entre premiers voisins dans le solide 
pur lorsque I=J et la moyenne entre celle de chaque constituant lorsque I 7^ J), et 
le rayon de coupure des interactions dont le choix est plus delicat. 
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3 Modele TB-SMA pour I'etude du CoPt 



Les interactions ne sont pas coupees brutalement, mais sont raccordees a zero a I'aide 
d'un polynome du cinquieme degre. Le raccord se fait ici entre la plus grande distance 
seconds voisins (Pt) et la plus petite distance troisiemes du (Co) de maniere a prendre 
en compte tons les premiers et seconds voisins pour chaque element dans I'expression 
de lenergie par site et pouvoir aussi negliger I'influence des troisiemes voisins de fagon 
continue et derivable. 

3.4.1 Energie de solution 

D'un point de vue theorique, I'energie de solution d'une impurete de A dans une 
matrice de B est definie comme la difference d'energie entre un etat initial constitue 
de Ns atomes de B dans un volume de B et de 1 atome de A dans un volume de A, et 
d'un etat final constitue de 1 atome de A dans une solution solide de B comportant 
Nfi-l atomes de B et de 1 atome de B dans un volume de B |6j. Cette energie 
s'exprime en utilisant I'expression de I'energie par site (2.2) de la fagon suivante oii 
les indices 1,2 et 3 representent respectivement les interactions de type Pt-Pt, Co-Co 
et Co-Pt: 



EcoiPt) = [-y^{l2Gq3du+QGq3d22)+l2Gp3d22+6Gp3d22]+e[-V^{l2Gqidii+5Gqid22+Gq3d22)+l2GpiC 

(3.5) 

avec: 

GqKdjj = eie-'''-^^-') 

Gpxdij = axe '^'''^ok 
avec: K=l,2,3. 

On peut exprimer de la meme fagon I'energie de solution d'un atome Pt dans une 
matrice de (Co). 

3.4.2 Energie de formation des phases LIq, A1 et LI2 de Talliage CoPt 

A partir de I'expression de I'energie d'une configuration (2.4), nous pouvons exprimer 
en remplagant les (Pf ) par leurs valeurs I'energie de formation de la phase LIq 
ordonnee de I'alliage CoPt: 



EfiLlo) = l/2*[{E\Pt)+E\Co)+E'{Pt)+E^{Co))-{E,oh{Co)+E,oh{Pt))] (3.6) 

avec I'energie d'attraction (terme de bande): 

E\Pt) = [-^{AGqidii + 8Gq3di2 + 4Ggid2l + 2Gqid22)] 

E\Go) = [-^{AGq2dii + SGq3di2 + 40^2^21 + 20^2(^22)] 

I'energies de repulsion entre les atomes: 

E'-(Pt)=(4Gpidii+8Gp3di2+4Gpid2i+2Gpid22) 
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3.5 Parametrage du potentiel.l 



W (Co)= (4Gp2dii+8Gp3di2+4Gp2d2i +2Gp2d22 ) 

et de la meme maniere, I'energie de formation de la phase Ai (LIq desordonnee) : 
Ef{Llo) = l/2*[{E\Pt)+E\Co)+E'{Pt)+E^Co))-{Ecoh{Co)+Ecoh{Pm (3.7) 
avec I'energie d'attraction: 

E\Pt) = [-^/l2Gqidu + 2Gq3du + "^Gqidu + ^Gqsdu + 2Gqid2i + 2Gg3d2i + 

Gqid22 + Gqzd22)] 

E\Co) = [-^/{2Gq2dn + 2Gq^dxx + AGqsdu + 4^92^12 + 2Gg2C?2i + 2^53(^21 + 
Gq2d22 + Gq3d22)] 

I'energies de repulsion entre les atomes: 

E^(Pt)=l/2*(2Gpidii+2Gp3dii+4Gpidi2+4Gp3di2+2Gpid2i+2Gp3d2i+Gpid2i+Gp3d2i) 
E'-(Co)=l/2*(2Gp2dii+2Gp3dii+4Gp2di2+4Gp3di2+2Gp2d2i+2Gp3d2i+Gp2d22+Gp3d22) 

De la meme maniere I'energie de formation de la phase LI2 ordonnee s'obtiennent: 

(LI2 ) = 1/4 * [E^i (Pt) + E^^ (Co) + 2 * (p^) ^ ^ri (p^) ^ ^r2 ^^^^ + 2*E'-^ {Pt)] - 

1/2 * [E^hiCo) + E^hiPt)] (3.8) 

Ou: 

I'energie d'attraction (terme dc bande): 
E^\Pt) = [-^{AGq^dii + 8Gqidi2 + ^Gqid2i + 2Gqid22)] 
E^^iCo) = [-^/lAGqadii + 8^93^12 + 4Gg2c/2i + 2092^22)] 
E^^Pt) = [-7(4^93^12 + 4Ggidii + 4Gqidi2 + 4Gqid2i + 2Gqid22)] 
I'energies de repulsion entre les atomes: 
E'-i(Pt)=(4Gp3dii+8Gpidi2+4Gpid2i+2Gpid22) 
(Co) =(4Gp3dii+8Gp3di2 +4Gp2d2i +2Gp2d22) 
E^3(pt)=(4Gp3di2+4Gpidii+4Gpidi2+4Gpid2i+2Gpid22) 

3.5 Parametrage du potentiel.l 

Nous decrivons dans cette partie la parametrisation des interactions atomiques croisees 
que nous avons effectue sur I'alliage CoPt les parametres (^3, A3, p3 et q3, I^^J) 

sont determines par I'ajustement des energies de formation des differentes phases du 
systeme CoPt ainsi que les distances interatomiques a I'equilibre de ces memes phases 
sur les valeurs experimentales. 

En faisant varier les parametres {^ij, A/j, p/j et qjj) dans plusieurs intervalles, nous 
avons calcule le parametre de reseau en minimisant I'energie de formation et nous 
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3 Modele TB-SMA pour I'etude du CoPt 



Par. 


Pot. 2 


Pot. 1 


A/j 


0.245 eV 


0.175 eV 


CiJ 


2.386 eV 


2.115 eV 


PiJ 


9.97 


9.412 


(liJ 


3.32 


2.812 



Table 3.2: Parameters du potentiel SMA pour I'alliage CoPt ajustes sur les valeurs 
experiment ales d'energie de dissolutions Pot. 2 |22] et de formation Pot.l 
fdes differents phases (LIq, Al et LI2) de I'alliage CoPt. 

avons determine cette energie de formation minimum qui a ete comparee a I'energie 
de formation experimental. 

Cette procedure est tres laborieuse, vu le nombre de parametres en jeu, mais en 
utilisant une approche numerique, nous avons pu par interpolation determiner un jeu 
de parametre reproduisant de fagon satisfaisante les valeurs experimentales. 

Les parametres retenus lors de la procedure d'ajustement numerique sont resumes 
dans le tableau 2.2, ainsi que ceux de Goyhenex et al. (Pot. 2). 

Nous avons illustre sur les figures(2.6, 2.7, 2.8) les resultats de I'ajustement des 
parametres croises, c'est-a-dire, revolution de I'energie de formation en fonction du 
parametre de reseau pour les deux potentiels, ceci dans le but de verifier que nous 
stabilisons bien les structures LIq, LI2 et Al pour le parametre de maille a egale au 
parametre experimental, et que d'autre part, notre potentiels equilibres ^ 

Les figures(2.6 et 2.7) montrent que notre potentiel Pot.l decrit bien I'energie de 
formation de LIq et Al, par contre le Pot. 2 ajuste un peu mieux le parametre de 
reseau experimental. 

L'accord est moins bon sur ce qui concerne la structure LI2 (figure 2.8). Neaumoins 
cette structure n'est pas I'objet de notre etude. 

L'ajustement des energies (tableau 2.4) sur les valeurs experimentales n'est pas 
parfait, mais c'est le meilleur que Ton puisse faire pour satisfaire toutes les conditions 
requises. 

Dans le tableau (2.3) sont rassemblees les valeurs experimentales des energies de 
formation des phases (LIq, Al et LI2), ainsi que les energies de solutions calculees a 
I'aide des deux potentiels, et les energies de solutions d'une impurete de Co dans une 
matrice (Pt) et Pt dans (Co), sur les quelles ont ete ajustes, les parametres (A/j, 
P/J) q/J) ?/j) des interactions Co-Pt. 

^partie attractive et repulsive d'egale importance 
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3.5 Parametrage du potentiel.l 



E (eV) 


Pot. 2 


Pot.l 


Exp 


E/Lio 


-0.2 eV 


-0.147 eV 


-0.145 eV 


Xi,io(% of a, c) 


1.002 


0.999 




Eai 


-0.125 eV 


-0.055 eV 


-0.10 eV 


Xai(% of a) 


1.007 


1 




E/L12 


-0.47 eV 


-0.41 eV 


-0.12 eV 


Xli2(% of a) 


1 


1.001 




EsoiiCo in Pt) 


-0.47 eV 


-0.33 eV 


-0.47 eV 


Esoi (Pt in Co) 


-0.64 eV 


-0.55 eV 


-0.64 eV 



Table 3.3: Energies de formation des differentes phases de I'alliage Co^ Pti-x et 
de solution d'un atome Co dans (Pt) et Pt dans (Co), et les valeurs 
experimentales equivalentes|22). 
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3 Modele TB-SMA pour I'etude du CoPt 



Type 


Eco«.(eV) 


E f-Ey 

J * 


Eco«.(eV) 


E f-Ej 


Saut d'atome 


sans relax. 




avec relax. 




Co — > Co 


1.74 





1.04 





Co — > Pt 


1.15 


0.245 


0.80 


0.10 


Pt — > Co 


1.735 





0.91 





Pt — > Pt 


0.58 


0.12 


0.30 


0.15 



Table 3.4: Energies de col calculees avec le Potentiel 1 

3.6 Calculs de barrieres d'energie (energie de col) 

Les barrieres d 'energies a franchir par un systeme lors de processus comme un saut 
atomique ou un glissement de dislocation, sont des grandeurs utiles a connaitre pour 
modeliser la cinetique d'un systeme, et en particulier lorsque I'on veut les simuler par 
une methode Monte-Carlo. La dtermination de I'energie de col dans un espace de 3N 
dimensions o N est le nombre d'atomes pouvant relaxer, est un problme difficile. Pour 
cela nous avons utilise la methode de Dynamique Moleculaire trempee en Liaisons 
Fortes avec I'approximation au second moment de la densite d'etats (TB-QMD) qui 
permet de chercher le chemin de plus grande pente reliant deux minima donns. Nous 
avons modifi I'algorithme de maniere a choisir un atome et de positionner une lacune 
premier voisin de cette atome, la methode utilisee est la suivante: 
On fait bouger cette atome vers la lacune et on suit revolution de I'energie totale du 
systeme en fonction de la position de I'atome qui saute, avec et sans relaxation (sans 
relaxation, on prend en consideration I'energie initiale du systeme). L'etude a etait 
faite sur les quatres types de sauts possibles, figure (2.9): 

-atome Co vers lacune sur site Co, 

-atome Co vers lacune sur site Pt, 

-atome Pt vers lacune sur site Co, 

-atome Pt vers lacune sur site Pt. 
A chaque pas de cet atome vers la lacune on calcule I'energie totale du systeme. 

Resultats 

Les courbes des figures(2.10 et 2.11) reprsentent I'energie totale du systeme CoPt 
de structure LIq calculee avec le potentiels interatomiques SMA Pot.l et Pot.2|^ 
en fonction de la position de I'atome qui saute, pour les 4 differents. 

L'energie totale du systeme augmente en fonction du deplacement de I'atome qui 
saute, puis devient maximale au milieu du deplacement, ensuite elle diminue lorsque 
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3.7 Conclusion 



Type 


Ecoi.(eV) 


E f-Ej 

J 


EcoZ.(eV) 


E f-Ej 


Saut d'atome 


sans relax. 




avee relax. 




Co — >Co 


2.48 


eV 


1.47 





Co — >Ft 


1.539 


0.318 


0.87 


0.27 


Pt — >Co 


2.12 





0.97 





Pt — >Ft 


0.66 eV 


0.144 


0.33 


0.13 



Table 3.5: Energies de col calculees avec le potentiel 2 
1' atome se rapproche de la position finale (lacune). 

On remarque que pour chaque saut la barriere d'energie est plus importante dans le 
cas sans relaxation, on peut done deduire que si on ne prend pas en compte les effets 
de relaxation on surestime les barrires de migration. 

On remarque aussi une legere difference entre I'etat initial et I'etat final dans le cas 
atome Co vers lacune Pt ou atome Pt vers lacune Co cette difference est due a la 
nature des fenetres^ de saut. Les resultats numeriques de ces energies de col sont 
rassembles dans les tableaux (2.5 et 2.6). 

3.7 Conclusion 

Dans ce chapitre, nous avons elabore notre modele energetique (Potentiel. 1) par 
ajustement aux donnees experimentales de I'alliage CoPt. Nous avons aussi presente 
le potentiel de Goyhenex et a/[^, (Potentiel.2). 

Nous avons applique ces potentiels pour le calcul des barrieres energetiques; necessaires 
pour une etude Monte-Carlo ulterieure dans les chapitre 3 et 4, ces potentiels seront 
utilises pour I'etude Monte-Carlo (MC-SMA) et en Dynamique Monleculaire (MD- 
SMA). 



^dans le cas d'un saut Co vers Co lacune, I'atome traverse une fenetre qui comporte 4 atomes de 
meme type Pt dans le cas d'un saut Co vers Pt lacune, I'atome traverse une fenetre qui comporte 
2 atomes de Pt et 2 atomes de Co 
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3 Modele TB-SMA pour I'etude du CoPt 




Figure 3.4: Maille representative d'une impurete de Co dans une matrice de (Pt) 
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3.7 Conclusion 
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3 Models TB-SMA pour Petude du CoPt 




-0.35 I ' ' ' ' ' ' ' ' ' 1 

0.9 0.95 1 1.05 1.1 1.15 

Variation des parameti'es de maille a et c (%) 

Figure 3.6: Variation de I'energie de formation de la phase LIq en fonction des 
parametre (a et c) de reseau 
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3.7 Conclusion 
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3 Models TB-SMA pour Petude du CoPt 



Efu.(eV) 



Systeme CoPt, (structure LL). 




Variation de parametre de maille a , (%). 



Figure 3.8: Energie de formation de la phase LI2 de I'alliage CoPt en fonction du 
parametre a de reseau avec les deux potentiel Pot. 1 et Pot. 2 
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3.7 Conclusion 




3 Models TB-SMA pour Petude du CoPt 
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A Sans relaxation 
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Figure 3.10: L'energie totale (eV) du systeme en fonction de la position de ratome 
qui saute vers la lacune (V), avec le Potentiel 2. 
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3.7 Conclusion 



tot 



A Sans relaxation 
-2635.0 - OAvec relaxation 
A 

A A 

-2636.0 

A A 

oOo 

-2637.0 - A A 
Co_Co(V) 

-2638.0 I ' ' ' 

-0.1 0.2 0.5 0.8 

Position relative 

-2605.0 



-2606.0 



-2634.5 



-2635.0 



-2635.5 



OA 13 
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-2608.0 





A 
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A A 


A 


oO°OoA 


A 

■ ^00 

1 
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Position relative 



-2636.0 F I Co_Pt(V) ""OA 

2636.5 
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-2631.7 
-2631.9 
-2632.1 F 

-2632.3 
-2632.5 
-2632.7 



^0 





Pt_Co(V) 



1.1 



-0.1 0.2 0.5 0.1 
Position relative 



1.1 



Figure 3.11: L'energie totale (eV) du systeme en fonction de la position de I'atome 
qui saute vers la lacune (V), avec le Potentiel. 1. 
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3 Modele TB-SMA pour Petude du CoPt 
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4 Transition Ordre-Desordre dans Talliage CoPt 



4.1 Introduction: 



De maniere generale, la simulation permet d'etudier un systeme donne dont on con- 
nait les interactions complexes, de mesurer les effets de certains cliangements dans les 
interactions sur le comportement du systeme et d'experimenter de nouvelle situations. 

Dans ce chapitre, nous nous interessons a la transition ordre-desordre de la phase 
equiatomique LIq de I'alliage CoPt. Pour cela on a repris le modele Monte-Carlo 
utilise auparavant avec un Hamiltonien de type Ising par Yaldram jHEl El dans la 
phase B2, E. Kentzinger en collaboration avec I'equipe de M. Benakki (LPCQ) j37ll38j 
dans les phases B2 et DOS, et repris par Oramus et al. [32], Kerrache et al.[10\ 
HO] dans les structures ordonnees sur les reseaux cubiques a faces centrees et M. 
Hamidi dans les alliages binaires de structure hexagonale compacte que nous 
modifions en introduisant les potentiels interatomiques calcules par I'approche plus 
realiste basee sur les Liaisons Fortes avec une approximation au second moment de la 
densite d'etats que nous avons decrite au chapitre precedent (Pot.l et Pot.2[?]). Nous 
commengons dans une premiere partie a presenter les differents algorithmes utilises 
dans les simulations. Nous decrivons ensuite notre modele theorique (MC-SMA) et 
nous presenterons les differents resultats obtenus dans I'alliage CoPt de structure 
Llo. 
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4 Transition Ordre-Desordre dans I'alliage CoPt 



4.2 Simulation Monte- Carlo 

L'etude d'un systeme comportant un tres grand nombre de degres de liberte se fait 
generalement en considerant un petit nombre d'atomes plonges dans un environ- 
nement moyen. Ainsi, les calculs des cinetiques de reaction en physique reposent sur 
I'approximation que chaque atome est entoure d'un nombre moyen de chacun des 
elements entrant dans la solution. L'utilisation de la methode Monte Carlo dans la 
simulation des transformations de phase a ete suggeree par Murray et utilisee pour 
la premiere fois par Fosdick |171 143 j . ces simulations on debute par des mecanismes 
d'echange direct entre les atomes. Flin et Mc Manus '44' furent les precurseur dans 
l'utilisation du mecanisme lacunaire dans les simulations des cinetiques de transfor- 
mation ordre-desordre dans des sructures cubiques. 

Dans notre etude, il s'agit de modeliser la diffusion atomique dans I'alliage CoPt 
ou les cinetiques sont decrites par un mecanisme lacunaire faisant intervenir des sauts 
d'atomes vers les lacunes EHl El EHl IIHI EHI , ce mecanisme etant reconu comme 
le mecanisme de base de la mobilite atomique dans les alliages binaires. 

La probabilite d'effectuer le saut d'un atome donne d'un site initial i vers un site 
final f depend de I'ecard energetique Ej-Ej: 



Les differents algorithme different dans la maniere de considerer la variations 
d'energie Ej — Ei lors du saut de I'atome. 

4.2.1 Algorithmes utilises dans les simulations 
4.2.1.1 Algorithme de Metropolis 

Get algorithme a ete I'un des premiers utilises dans les simulations Monte-Carlo en 
physique de la matire condense par Metropolis, en 1953 |21J. II est bien adapte au 
mecanisme d'echange direct dans la diffusions atomique bien qu'il soit utilisable dans 
le cas du mecanisme lacunaire et a des temperatures relativement elevees. 
Les principales etapes de I'algorithme avec un mecanisme d'echange sont: 

• Une lacune est choisie au hasard dans la boite de simulation et un atome proche 
voisins dans lequel sautera la lacune. 

• On calcule la variation d'energie AE associee au saut de I'atome dans la lacune. 
La probabilite de saut dans cet algorithme est donnee par la relation: 




(4.1) 



P(AE)=1 si A£; < 
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4.2 Simulation Monte-Carlo 



P(AE)=exp(^) si AS > 

• On precede au tirage d'un nombre aleatoire R uniformement distribue entre 
et 1. Si R> P, le saut a lieu, sinon le saut n'a pas lieu. Dans chacune de ces deux 
possibilites, la configuration d'atomes obtenue est consideree comme une nouvelle 
configuration. Elle est souvegardee pour le calcul des grandeurs thermodynamiques. 
Le temps physique peutetre associe au nombre de macro-pas ("Mcs"). 

L'algorithme de Metropolis permet tous les sauts d'atomes d'un etat initial vers 
un etat final d'energie inferieure (A E j 0), meme si entre les deux etats existe une 
barriere energetique. Pour A E ^ 0, revolution de systeme se fait avec une probabilite 
inferieure a 1. Get algorithme trouve sa plein efficacite a haute temperature. 



4.2.1.2 Algorithme de Glauber 

L'algorithme de Glauber suit les meme etapes que l'algorithme de Metropolis , 
mais le mode de calcul de la probabilite de saut les differencie. Dans cet algorithme 
la probabilite est donnee par la relation: 

Pi^E) = , '^"^^L (4-2) 

La probabilite d'effectuer des sauts pour AE < dans cet algorithme est toujours 
inferieur a 1, ce qui n'est pas le cas pour l'algorithme de Metropolis pT|. 

Pour AE =0, ou pour des temperatures elevees, le saut qui etait toujours accepte 
par l'algorithme de Metropolis, cette fois ci a autant de chance d'avoir lieu. 

On pent montrer que dans le cas de l'algorithme de Glauber ou pour l'algorithme de 
Metropolis, le systeme atteint un equilibre thermodynamique apres un grand nombre 
d'essais de sauts. 



4.2.1.3 Algorithme a temps de residence 

Le mode de calcul de la probabilte de transition est le meme que celui utilise dans 
les algorithme de Metropolis [22122] et de Glauber [^. La difference reside dans le 
mode du choix de I'atome qui fait le saut. Au lieu de choisir aleatoirement I'atome 
qui sautera vers la lacune, on calcule le nombre n d'atomes premiers voisins avec 
les probabilites correspondantes a chaque saut Wi, W2,. . . (de type Metropolis ou 
Glauber) . 

• On construit un segments de longueur L=^^^^ nWi comprenant n segments de 
longueur Wj chacun. 

• On tire au hasard un nombre aleatoire R entre et 1. 

• On choisit le saut vers le j atome si le produit R*L tombe dans le j '^'^^ intervalle 
de L en enfin, 
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4 Transition Ordre-Desordre dans I'alliage CoPt 



• on calcule le temps associe au saut t=-^. 

Dans cet algorithme, un saut a eu lieu pendant chaque cycle Monte-Carlo. Con- 
cernant le temps total, on peut soit ajouter a chaque pas le temps calcule (r), soit 
ajouter le temps t=r Ln(R') ou est un nombre tire au hasard; pour un systeme 
suffisamment grand les deux methodes sont equivalents [HSl EH ■ 

4.3 Simulation Monte-Carlo dans Talliage CoPt 
4.3.1 Description du modele (MC-SMA) 

Le modele est base sur un mecanisme d'echange entre la lacune et ses proches voisins. 
L 'etude a ete faite sur un reseau rigide de 32x32x32 atomes et une lacune choisie au 
hasard dans le cristal, en utilisant les conditions aux limites periodiques dans toutes 
les directions. 

Procedure utilisee 

Une lacune choisie au hasard dans le cristal, puis I'un de ses sites plus proches 
voisins est choisi au hasard. Si ce site est occupe par une autre lacune, on choisit un 
autre voisin. S'il est occupe par un atome, on calcule la variation d'energie AE de 
I'amas, associee au saut de cet atome de sa position initiale vers sa position finale. 

La probabilite de saut de I'atome dans la lacune est donnee par la probabilite de 
Glauber [SH: 

avec e^: energie de col correspend au saut, que nous avons neglige dans cette etude. 

La variation d'energie est calculee en sommant sur toutes les energies des atomes 
autour de la lacune et de I'atome qui saute avec un rayon egal a la distance en- 
tre second voisins, en faisant attention de ne pas compter deux fois les atomes de 
I'intersection de I'amas autour de la lacune et de I'amas autour de I'atome voisin 
choisi. 

La structure LIq de I'alliage CoPt choisie dans nos simulation est une structure 
cfc avec un parametre de maille moyen ac/c= 3.77 A. On a pas tenu compte de la 
tetragonalite de la structure LIq de I'alliage CoPt a cause du probleme lie a la varia- 
tion des parametres de maille (a et c) en fonction de la temperature figure (2.3, b) j24j . 
La transition ordre-desordre dans la structure LIq de I'alliage CoPt s'accompagne 
d'une transition structurale oii changement de symetrie cristalline de tetragonale a 
faces centres (tfc) a cubique a faces centrees (cfc) a la meme temperature. 
Ce parametre de maille moyen a^/c a ete deduit des courbes de revolution des energies 
des deux configurations des phases LIq et Al en fonction de leurs parametre de maille 
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4.3 Simulation Monte-Carlo dans Valliage CoPt 



Euo(eV/at) 



LL CoPt (Volume) 



-4.9 



-5 



-5.1 



-5.2 - 



-5.3 - 




-5.4 I ' ' • ' • ' • ' • ' ' ' ' 1 

3.5 3.6 3.7 3.8 3.9 4 4.1 4.2 

a^f^ (sans Tetragonalisation) 

Figure 4.1: L'evolution de I'energie de configuration de la phase LIO ordonne de 
I'alliage CoPt en volume en fonction du parametre de maille moyen &cfc 
sans tetragonalisation. 



ac/c (pris sans tetragonalisation). Les fig. (3.1, 3.2) montrent que I'energie de forma- 
tion des deux configurations LIq et Al est minimale a la distance ac/c=3.77 A, prie 
pour parametre de maille dans notre modele. 

4.3.2 Choix des parametres de simulation 

La dimension du systeme choisie est de 32 x 32 x 32 (16384 atome de Co et 16384 
atome de Pt), dans une phase equiatomique de structure LIq sans tetragonalisation, 
avec un parametre de maille moyens ac/c = 3.77A, le nombre de lacunes etant fixe 
a 1. Les energie de cols des atomes migrants ont ete prises egales a 0. les calculs 
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4 Transition Ordre-Desordre dans Valliage CoPt 



E^,(eV/at) CoPt (Volume) 

-4.74 



-4.84 - 



-4.94 



-5.04 - 



-5.14 



-5.24 - 



-5.34 

3.5 3.6 3.7 3.8 3.9 4 4.1 4.2 
a^.^^ (sans Tetragonalisation) 

Figure 4.2: L'evolution de I'energie de configuration de la phase LIO ordonne de 
I'alliage CoPt en volume en fonction du parametre de maille moyen ac/c 
sans tetragonalisation. 



ont ete fait par les deux series de parametres Pot. 1 et Pot. 2 [?] decrit au chapitre 
2, le potentiel 1 representant notre travail. L'energie du saut calculee en tenant 
compte de I'energie de touts les atonies voisins avec un rayon de coupure fixe aux 
deuxiemes voisins autour de la lacune et de I'atome choisie, en comptant les atomes 
de I'intersection des deux amas d'atomes une fois, avec une portee des interactions 
jusqu'aux second voisins. 

L'etat d'ordre de la structure est caraterisee a I'aide d'un parametre d'ordre a longe 
distance distance t/q^^)- 

Parametre OLD: 
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4.3 Simulation Monte-Carlo dans I'alliage CoPt 





Pot. 2 


Pot. 1 


exp. 


Tc 


580 K 


930 K 


1110 K 



Table 4.1: Resultats des simulations et experimentale de la transition ordre-desordre 
dans la phase equiatomique CoPt. 



'lOLD — ^ 

avec N nombre total d'atomes et et nombre d'atomes A sur les sous reseaux 
a et /3 respectivement. 

r\Qi^D= structure desordonnee 

^^OLD~ ^ structure completement ordonnee 



4.3.3 Resultats des simulations 

Nos simulations ont ete faites a la stoechiometrie AB. avec une lacune ce qui corre- 
spond a une concentration de lacunes C^=3. 05175 x 10"*^^, Pour une temperature 
donnee, nous suivons revolution de rjoiD avec le temp Monte Carlo le parametre a 
I'equilibre est note Vold^ I'etude est faite pour differentes temperatures, est revolution 
de rjQ^j-) nous permet de determiner la temperature critique pour les potentiels 
1 (notre etude) et 2. Le domaine est parcouru avec un pas de 10 K aux basses 
temperatures (domaines ordonnes) et un pas de 1 K a I'approche de la temperature 
critique. 

L'evolution du parametre d'ordre en fonction de la temperature figure (3.3) permet 
de deduire une temperature critique ordre desordre qui se situe a 930 K avec nos 
parametres et 580 K avec les parametres Pot. 2 [?] Tab. (3.1). La valeur obtenue par 
notre approche est voisine de Tc experimentale. 

Nous remarquons que tjqld reste eleve de OK jusqu'au voisinage de Tc puis chute 
brutalement a zero pour T = Tc, ce qui signifie que I'ordre dans les domaines ordones 
est eleve, ce qui montre que la transition est fortement du premier ordre. Ceci 
Concorde avec les resultats de Dahmani |24j . 
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4 Transition Ordre-Desordre dans Valliage CoPt 
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Figure 4.3: L'cvolution de paramctrc d'ordre a longue distance rj en fonction de la 
temperature, calcule par MC-SMA. 



4.4 Calcul de Tc a partir des energies de formation des 
phases LIq et Al du CoPt 

On peut envisagcr d'utiliser les resultats du modcle SMA (clip.2) tels que I'cncrgio 
de la phase ordonnee LIq et la phase desordonnee Al sur reseau cfc afin de determiner 
I'interaction effective de paire V entre premier voisins pour ensuite calculer la temperature 
critique correspondant a la mise en ordre de la phase LIq sur un reseau rigide. Ainsi, si 
on ne tient compte que des interactions effectives aux premiers voisins, I'energie de for- 
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4.4 Calcul de Tc a partir des energies de formation des phases LIq et Al du CoPt 



mation de la phase ordonnee LIq par rapport a la phase desordonnee Al s'exprime|ni: 

Elio - Eai = -V'f' (4.4) 

Les energies intervenant dans ces expression sont obligatoirement determinees sur 
reseau rigide. Ce reseau est choisi comme le reseau qui minimise I'energie de la 
phase ordonnee. Ainsi, les valeurs concernant les phases ordonnee et desordonnee de 
I'alliage CoPt sont notees dans le tableau 2.2 (chp.2). De cette etude on deduit que 
I'interaction effective de paire est egale a V'=^'== -0.07 eV avec le potentiel (Pot. 2) 
et -0.048 eV avec (Pot. 1). 

On pent comparer ces resultats a ce que donnerait la parametrisation de V deduite 
de I'energie de melange AHm definie comme la difference entre les energies de I'etat 
desordonne et celle de deux blocs de metaux purs dans les proportions de I'alliage. 
C'est ce qui avait ete fait dans une etude preliminairep^ en considerant les interac- 
tions aux premiers voisins et I'alliage cfc a I'equiconcentration c=0,5. L'energie de 
formation de ce melange qui vaut experimentalement -O.llOeV/at. a 1350 K j57j . 
s'exprime par: 

Ai7„ = -c(l - c)ZV^f'' = -SVf" (4.5) 

Ce qui conduit a une interaction de paires effective Y"^" de 0.037 eV/at. dont le 
signe positif indique bien une tendence a I'ordre mais qui est nettement plus faible 
en valeur absolue que les precedentes. 

On pent done en deduire les temperatures critiques des structures LIq et B2, 
obtenues en champ moyen pour une alternance de plans purs a partir de la formule: 

-2c(l-c)(Z-4Z0y 
Tc = — (4.6) 

oil Kb est la constante de Boltzman et Z' represente le nombre de liaisons inter- 
plans. Ceci conduit, dans le cas des structures LIq et B2, aux expressions suivantes: 

{W- pour la structure cfc 
^ , , 4.7 
-j^ pour la structure cc 

Les valeurs de la temperature critique de la phase equiatomque de I'alliage CoPt, 
deduite a partire des equations (3.7) sont sur le tableau (3.2). 

On remarque que nos parametres (Pot. 2) reproduisent bien la transition ordre- 
desordre (Llo=iAl) dans I'alliage CoPt en utilisant le modele[Hj base sur les energie 
de formation des deux phases ordonne et desordonne de I'alliage CoPt. 
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4 Transition Ordre-Desordre dans Valliage CoPt 





A Eform. (Pot.2) A Eform. (Pot.l) A Eform. (exp.) 


Tc exp 


Tc 


1740 K 1113 K 997.10 K 


1110 K 



Table 4.2: Temperature critique calculee par: Notre approche Pot. 1, Potential. 2 et 
les valeurs experimentales des energies de formation 



4.5 Conclusion 

Notre etude MC-TB-SMA, efFectuee avec le potentiel.l (notre travail), nous a permis 

dc calculcr la temperature critique de la transition LIq =i Al; temperature qui est 
en accord avec I'experience par rapport au potentiel 2, developpe en considerant les 
proprietes de I'alliage dilue (Coa;Pti_a;). 

Nous avons obtenu aussi un bon accord en utilisant une approche basee sur des 
potentiels de paires deduits des energies de formations. 
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5 Etude en surfaces de Talliage CoPt 



5.1 Introduction 

Le travail expose dans ce chapitre concerne les systemes bidimensionnels ou quasi 
bidimensionnels, c'est-a-dire les surfaces propres et les interfaces dans le domaine de 
la monocouche adsorbee. 

Sutface 

air 




Figure 5.1: Representation de la surface. 



Une large gamme de techniques experimentales [HHl (AES: Auger Electron Spec- 
troscopy [Sniinni) LEED: low Energy Electron Diffraction LEIS: Low Energy 
Ion Scattering 1)3, 64 , ISS: Ion Scattering Spectrometry, HEIS: High Energy Ion 
Scattering j65j . XPD: X-ray Photoelectron Diffraction 'HS', STM: Scanning Tunnel- 
ing Mivroscopy [HZI) FIM: Field Ion Spectroscopy|^£2j) XPS: X-ray Photoelevtron 
Spectroscopy [70], et plus recemment la diffraction de rayons X rasants[7J|) ont ete 
utilisees pour mettre en evidence et en particulier mesurer la segregation superficielle. 
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5 Etude en surfaces de Valliage CoPt 



On peut egalement se demander quel est le role des surfaces dans les transitions 
de phases volumiques (ordre-desordre ou demixtion-desordre)?, quel est le role de 
I'orientation cristallographique de la surface?, quelle est I'influence de la temperature 
et de la concentration volumique sur les phenomena?. A quelle profondeur s'etend la 
surface ?. 

Par exemplc Ic depot dc couches de Co ct Pt sur un substrat cntraine toujours un 
rearrangement de sorte que le Pt soit a la surface. Notre etude ayant pour but en 
particulier d'expliquer ce type de comportement. 

Dans ce chapitre, nous nous limiterons a quelques concepts fondamentaux. Apres 
une description des principaux modes de croissance possibles d'un adsorbat sur un 
substrat et les proprietes des couches minces d'alliages CoPt, nous exposerons notre 
etude qui porte essentiellement sur la modelisation par dynamique molculaire impli- 
quant un tres grand nombre d'atomes. Le potentiel atomique dans lequel evoluent 
les atomes est calcule par la methode TB-SMA exposee au chapitre 2. Nous etudions 
ensuite revolution de la distance (relaxation) entre le plan de surface (Pt) et de sub- 
surface (Co) de I'alliage CoPt suivant la direction [001], et nous evaluons I'energie 
d'adsorption d'un ou plusieur atomes de Co et Pt deposes sur I'alliage CoPt. En 
particulier, nous suivons revolution de cette energie d'adsorption en fonction de la 
dimension de marche deposee sur I'alliage CoPt. 

5.2 Proprietes structurales des surfaces et interfaces 

Definition 

Surface: Prontiere physique entre une phase liquide ou solide et une phase gazeuse 
ou un vide. 

Interface: Frontiere physique entre deux phases condensees (ex. solide-solide, 
solide-liquide, liquide-liquide) . 

Les zones superficielles d'un composant sont dircctement exposees aux agressions 
exterieures, et la surface est le siege de phenomenes specifiques: relaxation de la 
structure et reconstruction, presence de defauts, adsorption d'atomes etrangers . . . 
Elle peut differer considerablement de la structure du volume fig(4.2). Pour obtenir 
une surface, on coupe des liaisons interatomiques et les liaisons pendantes peuvent si 
licr entre elles et conduisent a une reconstruction de surface. On a alors une structure 
periodique avcc une nouvcllc periodicite; le dernier plan atomique peut simplement 
etre deplace par rapport a sa position normale dans le solide, c'est le phenomene de 
relaxation; Des atomes etrangers peuvent se fixer en surface, en des endroits precis 
(sites). 

Dans le cas d'un alliage, la concentration en surface de I'un des constituants peut 
etre superieure a cells du volume, c'est le phenomene de segregation. Dans un al- 
liage binaire A-rBi-j, on dit qu'il y a segregation a la surface de I'espece A si la 
concentration de A a la surface est superieure a x. Ce phenomene a des implications 
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5.2 Proprietes structurales des surfaces et interfaces 



00 00 00 : : : : § § § • • § 

00000 O^O^O^O o^o^o^o^o^o 

OOOOOOqoOO 

000000 • • • 

000000 0000 0000 000000 

Reconstruction Segregation Chimisorption Formation d'un composi 

Figure 5.2: Differentes configurations d'atomes au voisinage de la surface d'un cristal. 



importantes en metallurgie (corrosion, catalyse ...). La segregation de surface est 
fonction de la cristallographie (orientation de la face), la composition de I'alliage, 
la temperature, la presence d'impuretes residuelles en surface ou en volume. Les 
theories de la segregation ont mis en evidence deux facteurs principaux intervenant 
dans ce phenomene: 

- un efFet electronique du a la presence de liaisons coupees a la surface. 

- un effet elastique dii aux differences de rayons atomiques entre les deux constituants. 

Relaxation de surface 
Un atome de la surface n'a de voisin que dans son propre plan et en dessous. Com- 
pare a un atome du volume, il lui manque une force de maintien suivant la normale. 
On pent done prevoir que les plans atomiques superficiels sont deplaces suivant cette 
direction et II y a une modification de la structure des premieres distances inter- 
atomiques entre plans, la (fig.4.2) illustre cette situation. Toute fois I'amplitude de 
ces deplacements, souvent meme le sens, sont difficilement previsibles. 

Pour I'etude des reconstructions de surface ou des couches adsorbees qui forment 
un reseau periodique, on pent utiliser la diffraction d'electrons lents, la microscopie 
ionique de champ ou la microscopie tunnel. 



5.2.1 Modes de croissance 

La prediction theorique du mode de croissance d'un materiaux A (adsorbat) sur un 
substrat monocristallin S est un objectif delicat au regard du nombre de parametre 
intervenant dans les mecanismes de croissance cristalline. Les principaux modes de 
croissance possibles d'un adsorbat sur un substrat fig(4.4) sont: 

-Le mode de Volmer- Weber (VW): I'adsorbat ne mouille pas le substrat; on assiste 
a la formation d'ilots tridimensionnels des les premiers stades de la croissances (fig. 
4.4a); un exemple type de ce mode est donne par le systeme Pb/graphite [72j . 
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1 























• 


• 


• 


• 


# 




2 































-8- 








3 


-i- 


-§- 




-§— 








4 




• 


• 


• 


• 


• 


• 


5 



Figure 5.3: Representation shematique d'une relaxation de la surface: en blan,les po- 
sitions des atomes avant relaxation; en noir, les positions apres relaxation. 



-Le mode Stranski-Krastanov (SK): I'adsorbat mouille le substrat sur une ou 
plusieurs couches monoatomiques, et au-dela d'un certain recouvrement, la crois- 
sance se poursuit sous forme d'ilots (fig.4.4b); exemple Pb/Ge. 

-Le mode Frank- Van der Merwe (FVdM) fig(4.4c): la croissance se fait par couches 
successives, chaque couche etant remplie avant le debut de la croissance de la suivante; 
c'est le cas notamment du systeme Co/Cu(100)|2Hl Ell- 

Pour elaborer ces depot il existe plusieurs techniques, on distingue deux grandes 
categories: 

(i) - les methodes physiques telle que la pulvarisation ou I'evaporation, 

(ii) - les methodes chimiques telle que la CVD (chemical vapor deposition) et le 
depot electrochimique. 

5.2.2 Les defauts de surface 

II n'y a pas de raison pour que la plupart des defauts existant dans le volume n'aient 
pas leurs equivalents a la surface. On retrouve effectivement lacunes, interstitiels, 
impuretes ainsi que leurs associations, emergences de dislocations .... II y a tout 
lieu de penser que ces defauts jouent un role dans les phenomenes d'adsorption, de 
nucleation, qui sont determinants dans la croissance par epitaxie. 
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5.3 Proprietes des surfaces de Valliage CoPt 
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(a) (b) ■ (c) 



A 7 7 7 V 7 r 7 7 r ~ ~ r 7 7 ^ 7 7 r 



Volmer Weber Stranski-Krastanov Frank-van derMerwe 

Figure 5.4: les trois principaux modes de croissance 

5.3 Proprietes des surfaces de Talliage CoPt 

De nombreuses etudes effectuees sur les multicouches out mis en evidence le role 
determinant des interfaces Co/Pt dans I'obtention de materiaux avec une anisotropie 
magnetique perpendiculaire [T^ • H n'est pas question ici de faire une bibliographie 
extensive sur ce vaste sujet. L'article de revue d'Yves Gautier [JHI presente un bilan 
des etudes de segregation dans plusieurs series d'alliages monocristallins a base de 
Pt. En resume, on pent dire qu'apres quelques annees consacrees aux etudes sur les 
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5 Etude en surfaces de Valliage CoPt 



I 1 



II 



SM 



MSM 



Figure 5.5: croissance en multicouches simultanees (SM), croissance monocouche puis 
multicouches simultanees (MSM) 



multicouches, le role tres important de I'anisotropie magnetique des interfaces Co- 
Pt pour le developpement d'une anisotropic perpendiculaire et d'un champ coercitif 
eleve dans ces mater iaux en couches minces |75| a ete mis en evidence. Les chercheurs 
se sont alors orientes vers la preparation de la multicouche ideale: alternance d'une 
monocouche de Co et d'une monocouche de Pt, qui se trouve etre la structure de la 
phase Llo le long de sont axe tetragonal (001). 

Les premiers travaux sur les materiaux " artificiels''' CoPt sont, ceux de Lairson 
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5.4 Relaxation de surface dans le CoPt 



et al\n\ [THl ES]. lis ont etudie les multicouches Co/Pt preparees par pulvarisa- 
tion cathodique a 370c sur substrat de Mgo (001) dans le but d'obtenir la phase 
LIq ordonnee. La meme equipe a prepare ulterieurement des alliages CoPt par co- 
pulverisation a difFerentes temperatures. lis ont constate que le parametre d'ordre 
augmente avec la temperature de preparation. 

La technique de preparation par epitaxie par jet moleculaire a aussi ete utilisee avec 
succes pour co-deposer des couches minces d'alliage CoPt. Harp et ces collabor ateurs [50] 
ont depose plusieurs series d'echantillons de lOOOA d'epaisseur a temperature vari- 
able sur une couche tampon de Pt. Avant de passer a I'etude des depots d'atomes 
de Co et Pt sur I'alliage CoPt nous presenons dans le paragraphe suivant la distance 
entre la plan de surface (Pt) et de subsurface (Co) apres relaxation avec le modele 
dynamique moleculaire base sur la Liaisons Fortes. 



5.4 Relaxation de surface dans le CoPt 

Pour evaluer revolution de la distance entre le plan de surface et de subsurface apres 
relaxation dans I'alliage CoPt (cfc), on a utilise le modele de dynamique moleculaire 
presente dans le chapitre 2, avec comme potentiel d'interaction le potentiell. Les 
simulations on ete faites sur un alliage equiatomique LIq CoPt avec les conditions 
aux limites periodiques suivant les directions ox et oy et les conditions limites fixes 
suivant la direction oz. 

Cette distance entre le plan de surface et de subsurface est donnee par la moyenne 
des distances interatomiques suivant la direction [001], entre les atomes de surface 
(Pt) et subsurface (Co) (Fig. 4.7) s'ecrit: 

z=^^ (5.1) 

avec N et le nombre de paires Pt-Co, et Xj, la distance entre Pt-Co suivant la 
direction [001]. 

Pour le calcul de cette distance on a pris un reseau de structure (cfc), de 10 atomes 
dans les directions X et Y avec les conditions aux limites periodiques, et une dimen- 
sion variable suivant Z avec des conditions aux limites fixes. 



Cette distance avant relaxation vaut z=a/2 =1.885 A, (a parametre moyen de la 
phase cfc egal 3.77 A). On laisse ensuite le systeme se relaxe et on calcule a nouveau 
cette distance dans les cas: 

- d'un systeme avec 10 atomes suivant la direction [001], 

- avec un systeme de 20 atomes dans la direction [001], 

- et un systeme qui comporte 30 atomes dans cette direction. 

Ceci dans le but d'etudier I'influence de I'epaisseur sur la relaxation. Les resultats 
sont donnes dans le tableau (4.1). 
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5 Etude en surfaces de Valliage CoPt 




Figure 5.6: Boite de simulation utilise (avec (Pt) en surface) 



Considerons maintenant le cas de la distance entre un atome Pt du plan de surface 
avec son premier voisin Co de subsurface avant relaxation: x=a /y^2)=2.665 a. 
Apres relaxation, les valeurs trouvees pour les trois epaisseurs considerees sont donnees 



Table 5.1: Distances entre le plan de surface et de subsurface de I'alliages CoPt; apres 
relaxation. Sans relaxation z=1.885 A 



10 atomes 


20 atomes 


30 atomes 


1.7615 A 


1.7620 A 


1.7600 A 
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5.5 Adsorption 
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Figure 5.7: distance entre le plan de surface et de subsurface 



dans le tableau 4.2. 

On remarque que ces distances ne changent pas en fonction de I'epaisseur et que la 
relaxation induit une contraction de 6.63 % des distances entre les plans de surface 
et de subsurface. 

5.5 Adsorption 

L'adsorption est une technique efficace pour obtenir des informations sur la nature 

de la surface, pour estimer I'aire de surface et le degre d'adsorption, etc . . . 

Les molecules et les atomes peuvent se fixer sur une surface d'un substrat de deux 
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5 Etude en surfaces de Valliage CoPt 



Table 5.2: Distances entre un atome Pt du plan de surface avec son premier voisin 
Co de subsurface, apres relaxation. Sans relaxation x= 2.665 A 



10 atomes 


20 atomes 


30 atomes 


2.5799 A 


2.5800 A 


2.5800 A 



f aeons: 

• adsorption physique (pliysisorption) : I'interaction entre la molecule ou I'atome 
(adsorbe) et la surface (substrat) est faible; -AH jlOOkJ/mol. La variation d'enthalpie 
est done trop faible pour provoquer une rupture de liaison. 

• adsorption chimiquc (chimisorption): les molecules ou les atomes adherent a la 
surface en formant une liaison chimique (nature covalente) et le nombre de coordi- 
nation avec le substrat est maximum. La distance entre la surface et I'adsorbat est 
plus courte et on a - AH ^lOOkJ/mol. 

On pent aussi definir le degre d'adsorption d'une surface couverte d'atomes ad- 
sorbes par: 



nb. de sites d'adsorption occupes 
nb. de sites d'adsorption disponibles 

r est appele le taux de recouvrement et depend de la pression du gaz 
et la variation de r avec la pression est appelee I'isotherme d'adsorption. 



5.5.1 Energie d'adsorption d'un atome Pt ou Co 

Pour calculer le gain d'energie correspondant a I'adsorption d'un atome depose sur 
I'alliage CoPt, Tine etude par simulation Dynamique Moleculaire en tenant compte 
de la relaxation a ete entreprise. 

Pour calculer I'energie d'adsorption, il faut comparer I'energie totale des deux systemes 
representees sur la (Fig. 4.8) oii un atome a ete ajoute sur la surface. Cette energie 
d'adsorption s'ecrit: 

EadsA^ = Etot(A + b) - Etoi(A) - mu(Co) 

oii: 

- b atome depose sur A 
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5.5 Adsorption 



- A Matrice (alliage CoPt) 

- mu(Co) est le potentiel chimique de Co dans sa phase vapeur. On le prend en 
general comme origine des energies, egal 





Figure 5.8: 



Ces energies EadsA'' et Etot(A + b) ont ete calculees a I'aide de notre modele TB- 
SMA-QMD en utilisant notre potentiel (Pot.l). les valeurs trouvees sont donnees sur 
le tableau 4.3. 

II n'existe pas a notre connaissance de valeurs experimentales avec lesquelles on 
pourrait comparer nos resultats. Neanmoins le signe ainsi que I'ordre de grandeur des 
energies d'adsorptions calculees sont corrects sont du mme ordre de grandeur que les 
nergies de cohsion des mtaux purs (ce qui est quand mme bon signe) et qui plus est 
des valeurs comprises entre I'nergie de cohsion de Co et celle de Ft (Eco/i.(Ft)=-5.860 
ev et Eco/i.(Co)=-4.450 ev). 

On pent par contre citer les resultats des energies d'adsorption de C. Goyhenex el 
al. pour: 

Co/Ft(lll)[?]: E,rf,= -4.65 eV. 

Ft/Co(001)|ll|: E„rf,= -5.37 eV. 

Valeurs du meme ordre de grandeurs que les notres. 
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5 Etude en surfaces de Valliage CoPt 




5.5.2 simulation de marches de Co ou Pt sur I'alliage CoPt 

Nous considerons maintenant le depot de marches depaisseurs monoatomique sur 
I'alliage CoPt (cfc). Nous pouvons definir I'energie d'adsorption par atome de cette 
marche par: 

EtotiO = Nco,Pt) - Etotie = 0) , ^. 

Eads = Q (5.3) 

oil 9 est le nombre d'atome constituant la marche. 



Table 5.3: Les energies d'adsorptions d'un atome Pt et Co depose sur I'alliage CoPt, 
apres relaxation. 



Eafes(Pt) 


Ea6s(Co) 


-4.9951 eV/at 


-4.4649 eV/at 
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5.5 Adsorption 




Figure 5.10: Marche de 40 atomes de Pt sur I'alliage CoPt. 



Dans cette etude, nous avons simule a I'aide de notre modele (QMD-TB-SMA) et 
en utilisant notre potentiel (Pot.l), le depot de marche de Co ou Pt sur la surface 
(Pt) de I'alliage CoPt (cfc); Le nombre d'atomes constituant la marche variant en- 
tre 10 =i 100 atomes (100 atomes correspend a une couche complete) fig(4.10 et 4.11). 

Les resultats obtenus sont donnes dans le tableau 4.4 et dans la figure(4.12) ou 
nous avons aussi inclu I'energie d'adsorption d'un seul atome. 

Comme pour le cas de I'energie d'adsorption d'un atome, nous ne pourons pas com- 
parer nos valeurs a des resultats experimentaux. L'energie d'adsorption par atome 
diminue en valeur algebrique lorsque la dimension de la marche augmente dans les 2 
cas (marche de Co ou de Pt). Ceci indique que I'adsorption d'un atome est facilitee 
par la presence d'une marche de dimension importante. 

Cette energie est plus faible dans le cas d'une marche de platine; suggerons aussi 
que I'adsorption de Pt sur la surface Pt de I'alliage CoPt est plus aisee que celle du 
Co. 
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5 Etude en surfaces de Valliage CoPt 




Figure 5.11: Marche de 20 atomes de Co sur I'alliage CoPt. 



5.6 Conclusion 

Nous avons presente dans ce chapitre I'etude faite par dynamique moleculaire (QMD- 
TB-SMA) de surface de I'alliage CoPt. La relaxation de surface obtenue a ete evalue a 
6.6%. Ensuite, nous avons pu calculer les energies d'adsorption de differents adsorbats 
sur la surface Pt de I'alliage CoPt; le cas de la surface Co sera etudie ulterieurement. 

L'etude a ete faite uniquement en utilisant les parametres de notre potentiel 
(Pot.l). 
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5.6 Conclusion 



Nobre. d'atomc 


Nature de la marche 




de la marche( 9) 


Co 


Pt 


10 


-4.6024 


-5.4456 


20 


-4.7266 


-5.6633 


30 


-4.7489 


-5.7232 


40 


-4.7610 


-5.7535 


50 


-4.8281 


-5.7713 


60 


-4.7720 


-5.7820 


70 


-4.7760 


-5.790 


80 


-4.7780 


-5.7965 


90 


-4.7800 


-5.800 


100 


-4.7970 


-5.842 



Table 5.4: Energies d'adsorption par atome (en eV) de differentes marches de Pt et 
Co depose sur I'alliage CoPt. 
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Figure 5.12: 
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6 Simulation MC-lsing dans les couches minces 



6.1 Introduction 



Cette etude presentee a ete entreprise en parallele avec I'etude sur I'alliage CoPt 
(TB-SMA) dans le but d'apporter un complement au travail consequente effectue au 
sein de notre laboratoire et qui a porte sur les simulations des proprietes statiques 
(diagramme de phase) et cinetiques (energie de migration atomique) en utilisant un 
Hamiltonien de type Ising avec une portee des interactions aux seconds voisins. Ces 
travaux ont portes sur les transitions D03 — > B2 ^ A2 (cc) [?], LIq — > Ai (cfc) [?] 
et D019 Bi9 — > A3 (hep) [?]. Toutes ces etudes ont ete effectues avec une boite de 
simulation de taille variable et des conditions aux limites periodiques dans les trois 
directions pour simuler le volume. 

Une etude simulant une couche mince de structure LIq de largeur M par I'introduction 
des conditions aux limites fixes dans une direction a ete jugee interessante a entre 
prendre dans le but d'etudier I'influence de I'epaisseur sur les proprietes physiques 
de la couche mince; et en particulier verifier que I'on retrouve les valeurs du volume 
lorsque M est grand et que Ton obtient celle du reseau carre (2D) ^] lorsque M 
diminue. 
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6.2 Description du modele 

Dans cette etude, nous avons repris le modele utilise auparavant dans notre labo- 
ratoire, pour plus de details, se referer aux theses de A.Kerrache [?] et M.Hamidi 

[?]■ 

L'Hamiltonien utilise est de type Ising, et les potentiels d 'interactions de paires ont 
ete choisis tels que V^^ = V^^ = -V^^, avec une portee des interactions jusqu'aux 
deuxiemes voisins (i=l,2). 

Les limites de stabilite des phases ordonnees sont determinees par les interactions 
de paires effectives = l/A(V\^+ Y^b-'^^ab) 

La probabilite de saut de I'atome dans la lacune est donnee par la probabilite de 
Glauber [3]: 



energie de col. 

Notre choix s'est porte sur cette probabilite plutot que celle utilisee dans I'algorithme 
de Metropolis jHH pour qu'a hautes temperature (etat desordonne) ou pour AE = 0, 
le saut ait la meme probabilite d'avoir ou de ne pas avoir lieu; de plus aucun saut 
n'a une probabilite egale a 1 lorsque AE reste finie. 

Un nombre entre et 1 est tire au hasard. Si ce nombre est plus petit que P{AE) 
le saut a lieu, si non il n'a pas lieu. L'unite de temps utilise est le cycle Monte-Carlo 
qui correspond a un nombre d'essais de saut egal au nombre de lacunes. 



6.3 Choix des parametres de simulation 

Les simulation sont effectuees dans un reseau de structure LIq, de 64 cellules dans 
les directions X etY, avec des conditions aux limites periodiques, et une dimen- 
sion variable (8, 16, ....,192 cellules) suivant Z avec des conditions aux limites fixes. 
Chaque simulation est entreprise de 8 a 10 fois (selon qu'on s'approche ou non de la 
temperature critique Tc) dans le but d'ameliorer la statistique et de tester la repro- 
ductibilite des resultats. 

Les energies d'interactions de paires effectives entre atomes premiers voisins Vi a 
ete fixe a 0.02 eV et second voisins V2 variant entre (-0.002 a -0.01). Les energies 
de col des atomes migrants ont ete prises egales a zero (e^ = = 0) ainsi que les 
energies d'interactions atome-lacune (V Av = ^Bv =0). Une une seule lacune dis- 
tribuee initiallement au hasard, (C^=cst). 
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6.3 Choix des parametres de simulation 



Debut 



V 

Choix d'un site 
proche voisin 



M 



Le site est une lacune 



Calculer I'energie de 
transition 



o 



f 

Calculer la probabilite: P 



f 

Generer un nombre aleatoireR 




Accepter le saut 



Figure 6.1: Organigramme de I'algorithme de simulation 
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6 Simulation MC-Ising dans les couches minces 



6.4 Resultats 

6.4.1 Calcul des temperatures critiques 

Dans cette partie, nous determinerons les temperatures critiques associees aux tran- 
sitions Llo Ai pour differentes valeurs de V2 et d'epaisseur de film. 

L'etat d'ordre de la structure est caracterise, de la meme fagon qu'au chapitre 3, 
par un parametre d'ordre a longue distance: 

eq. _ 2fM__^ 
^OLD — ^ 

avec N nombre total d'atomes et et nombre d'atomes A sur les sous reseaux 
a et /3 respectivement. 

''10LD~ *-* structure desordonnee 

VoLD~ ^ structure completement ordonnee. 
Pour VI et V2 fixes, nous suivons revolution de rjQj^j-) avec le temps Monte Carlo, 
pour unc temperature donncc. 

Les courbes obtenues sont decrites par une loi de type exponentielle 

ri{t) = ai + a2exp{-t/rjL) 

L'ajustement de ces courbes avec cette loi nous permet de determiner le parametre 
d'ordre a I'cquilibre iJold ^ ™® temperature donnee ainsi que le temps dc relaxation 
t(T) correspondant. Un exemple est donne dans la figure5.2 pour VI = 0.002 eV et 
V2 = -0.004 eV. 

). En faisant varier la temperature, nous avons suivi revolution de rjoLD POur 
differentes epaisseurs est potentiels consideres et les resultats obtenus (figure5.3) nous 
permettent dc dcduirc les differentes Tc, ces valeurs seront discutccs par la suite. On 
remarque, a travers les allures des differentes courbes pour les differentes epaisseur 
que la nature de la transition ordre-desordre est du premier ordre. 



6.4.2 Energie de migration 

L 'evolution des temps de relaxation en fonction de la temperature est donne dans 
la courbe d'Arrhenius (fig5.4). Ces courbes sont ajustces a I'aide de la relation 
Ln(r) = Ln{To) + Em/KbT (avec Em est I'energie de migration). Ceci nous permet 
d'obtenir, pour chaque cas envisage, la valeur de I'energie de migration atomique Em 
les resultats sont donnes dans les tableaux 5.1 (M=48) 5.2 (M=64) 5.3 (M=128). 

On remarque que les parties lineaires des courbes pour de faibles epaisseurs sont peu 
etendues avec beaucoup de fluctuations comparativement au tres grandes epaisseurs, 
ce qui rend plus difficile la determination des energies de migration. 
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0.6 ' ' ^ ' ^ ' ^ ' ^ ' ' 

20000 40000 60000 80000 le+05 

Temps Monte-Carlo 

Figure 6.2: Variation du Temps Monte-Carlo en fonction de la temperature 



6.4.3 Discussion 

La figure 5.5 presente la correlation existant entre I'energie de migration et la temperature 
critique. Dans les etudes en volume [?], la dependance obtenue etait lineaire. Dans 
notre cas, on observe cette linearite pour M=64 et 128. Pour les faibles epaisseures, 
I'importance des fiuctuation est a soulignier. 

D'autre part, revolution de Tc en fonction de V2 est quasi fineaire pour toutes les 
epaisseurs (fig 5.6). Ceci reproduit les resultats en volume[?]. 



V2 


-0.002 


-0.004 


-0.008 


-0.01 


V2/V1 


-0.1 


-0.2 


-0.4 


-0.5 


Tc(K) 


537 


669 


932 


1052 


Em (eV) 


0.257 


0.303 


0.3279 


0.3455 



Table 6.1: Resultats de I'ajustement des courbes de Ln(rL)=f(l/T) pour I'epaisseur 
M=48 couches 
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Figure 6.3: Variation du parametre d'ordre a longue distance avec la temperature 
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Figure 6.4: Variation de Ln(rL) en fonction de 1000/T pour differentes epaisseurs. 
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V2 


-0.002 


-0.004 


-0.008 


-0.01 


V2/V1 


-0.1 


-0.2 


-0.4 


-0.5 


Tc(K) 


542 


676 


938 


1060 


Em (eV) 


0.258 


0.282 


0.339 


0.354 



Table 6.2: Resultats de Tajustement des courbes de Ln(rL)=f(l/T) par la methode 
des moindres carres avec une loi de type Ln(rL)= Ln(ro)-|- Em/K^T pour 
I'epaisseur M=64 couches 



V2 


-0.002 


-0.004 


-0.008 


-0.01 


V2/V1 


-0.1 


-0.2 


-0.4 


-0.5 


Tc(K) 


550 


682 


940 


1070 


Em (eV) 


0.2360 


0.2739 


0.3183 


0.3780 



Table 6.3: Resultats de I'ajustement des courbes de Ln(ri)=f(l/T) pour I'epaisseur 
M=128 couches 



• Influence de I'epaisseur du film: 

On a represente dans la figure 5.7 revolution de Tc en fonction de I'epaisseur pour 
les differents potentiels utilises. 

L'allure des differentes courbes est similaires, la temperature ordre-desordre augmente 
en fonction de I'epaisseur du film et a partir d'un certaine epaisseur elle presente 
clairement une valeur limite a partir de M=64 monocouches. Cette valeur limite 
coincide avec celle obtenue dans I'etude en volume de Kerrache [SJj pour une boite 
de L2=64x64x64 cellules avec des conditions aux limites periodiques, pour 

les meme valeurs de V2 (-0.004, -0.008, -0.01 eV). On pent dire done qu'a partir de 
cette epaisseur on obtient les proprietes du volume. 

Les energies de migration Em correspondentes aux differentes epaisseurs et po- 
tentiels sont rassembles dans la figure(5.8). D'apres cette figure on remarque que 
les traces d'Arrhenius dans le cas limite des faibles epaisseurs, nous donnent des 
energies de migrations qui coincident avec les resultats publies dans le cas du reseau 
bidimensionnels [HJ| . 



94 



6.4 Result ats 




Les resultats de Ejv/ dans le cas limite des grandes epaisseurs sont similaires avec 
obtenus dans le volume de structure LIq [?], figurc(2.15), ce qui veut dire qu' a partir 
d'un certaine epaisseur du film on attaint les proprietes du volume. 
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6 Simulation MC-Ising dans les couches minces 




400 ' • ' • ' • ' • ' • ' ' 

-0.6 -0.5 -0.4 -0.3 -0.2 -0.1 

Figure 6.6: Variation lineaire des temperatures critiques en fonction des energies de 
paires effectives 



6.5 Conclusion: 

Nous avons etudier dans le cadre de ce travail, par simulation Monte-Carlo I'interdependance 
de I'ordre et de la mobilite atomique dans les couches minces metalliques d'alliages 
binaires de structure LIq. 

Un mecanisme d'echange atome-lacune a ete utilise pour simuler les processus de 
diffusion dans I'alliage en ne prenant en compte que I'interaction atome-atome. 

Les cinetiques d'ordre a longue distance presentent un comportement de type ex- 
ponentiel. La transition ordre-desordre est plus ou moins du premier ordre. Nous 
avons trouve que la temperature critique augmente avec la taille du film. Elle se 
stabilise a partir de 64 couches atomiques ce qui determine la phase volume [HH]. Les 
temperatures de transition suivent une variation lineaire en fonction des energies de 
paires effectives pour chaque epaisseur du film. Les tracs d'Arrhenius montrent que 
les energies de migration atomique sont plus elevees aux faibles epaisseurs, et cela 
est du a I'effet de surface. 
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Figure 6.7: Variation de la temperature critique en fonction de la taille du film 
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Figure 6.8: Correlation entre I'energie de migration Em et I'epaisseur du film 
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6 Simulation MC-Ising dans les couches minces 
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7 Conclusion generale 



Dans ce memoire, nous avons presente I'elaboration d'un modele energetique decrivant 
les interactions interatomiques dans le cadre du formalisme des Liaisons Fortes ou 
la densite electronique est ecrite dans I'approximation du second moment. Pour 
decrire le systeme CoPt, nous avons ajuste les parametres du modele decrivant les 
interactions Co-Pt sur I'energie de formation des structures LIq et Al, ainsi que sur 
leurs parametres de reseau. Les parametres qui caracterisent les interactions Co-Co 
et Pt-Pt ont ete reprit de la litterature. 

Nous avons aussi utilise, dans un but de comparaison, le potentiel TB-SMA parametre 
par Goyhenex et al. [121 sur les energies de dissolution du Co dans (Pt) et du Pt dans 
(Co). Partant d'un code de dynamique moleculaire mis a notre disposition par le Prof. 
Treglia j84j de I'Universite de Marseille, nous avons evalue les energies de col associees 
aux differents sauts atomes lacunes; les valeurs obtenues montrent que la migration 
atomique est favorisee par la relaxation. D'autre part, nous avons estime les energies 
d'adsorption de differents depots de Co et Pt sur la surface de I'alliage CoPt ordonne 
avec comme derniere couche du Pt. Nos resultats montrent que I'adsorption de Pt 
est favorisee par rapport au Co et que I'energie d'adsorption diminue en fonction 
du nombre d'atomes adsorbes. A defaut de valeurs experimentales, nous pouvons 
uniquement dire que I'ordre de grandeurs des energies d'adsorption est correcte. 

En utilisant une approche Monte-Carlo, nous avons evalue la temperature critique 
de la transition LIq — > Al en utilisant notre potentiel ainsi que le potentiel de Goy- 
henex et al. j22j. La valeur obtenue avec notre potentiel est de 930 K, par contre en 
utilisant le second potentiel^^ nous obtenons Tc=580 K. Nous pouvons done con- 
clure que notre modele, parametre sur les proprietes de I'alliage equiatomique donne 
une meilleure description de Tc. Ceci a ete confirme dans une aproche utilisant des 
potentiels de paires deduits des energies de formation LIq et Al calcules a I'aide de 
notre potentiel. 

Dans ce memoire, nous avons egalement presente en parallele une etude utilisant 
un modele energetique different base sur des potentiels de paires effectifs, etude qui 
s'inscrit dans la continuite des simulations Monte Carlo realisees au sein de notre 
laboratoire. Ce travail a montre qu'en utilisant des conditions aux limites fixes dans 
une direction (dans le but de simuler une couche mince) on retrouve les resultats 
correspondant au volume lorsque I'epaisseur de la couche est grande {l 64 mono- 
couches) et de plus, lorsque I'epaisseur diminue (j Smonocouches) on tend vers les 
valeurs obtenues dans I'etude du reseau carre (2D) (HHI 

En perspectives, nous envisageons: 

- D'integrer les energies de col obtenues par dynamique moleculaire dans le code 
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7 Conclusion generale 



de simulation Monte Carlo dans le but d'obtenir les energies de migration atomique 
en volume. Soulignons que I'introduction d'energie de col ralenti considerablement 
les simulations Monte Carlo. 

- D'evaluer par dynamique moleculaire les energies de col associees aux divers 
sauts possibles en surface dans le but d'etudier par simulation MC la croissance et 
la morphologic des surfaces. 
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